= 
— 
7 
im 
: 
4 
~T 


44 
q 
4 
| 
» 
> 
i 
is 
i 
of Applied | 
“ 
‘ 
ag 
3 
= : 
= 
3 


- 
a 
- 
- = 
q 
g 
q 
q 
i 
q 
is 4 
j 
4 
4 
3 
: 
} 
+ 
: 
4 
q 
4 
- 
| 
| 
om 
a 
q 
"Be 


Veer 


OF 


EDITED BY 


H. L. DRYDEN T. C. FRY 
J. M. LESSELLS W. PRAGER 
J. L. SYNGE 


WITH THE COLLABORATION OF 


H. BATEMAN M. A. BIOT 

J. P. DEN HARTOG H. W. EMMONS 

J. N. GOODIER F. D. MURNAGHAN 
R. V. SOUTHWELL G. I. TAYLOR 


QUARTERLY 


APPLIED MATHEMATICS 


TH. v. KARMAN 
I, S. SOKOLNIKOFF 


L. N. BRILLOUIN 
K. O. FRIEDRICHS 
W. R. SEARS 

S. P. TIMOSHENKO 


VotumE 


1943 


a 
4 
| 
a 
| 
7 a 


Printed by the 
GEORGE BANTA PUBLISHING CoMPANY 
Menasha, Wisconsin 


: 
4 

. 

ap 

Bex. 


CONTENTS 


W.S. Ament: The lines of principle stress in the plane problem of plasticity . 

H. Bateman: The transformation of partial differential nennane ; 

L. Bers: Concerning the acceleration potential , 

and A. Gelbart: On a class of differential equations in , mechanics of 
continua 

L. Brillouin: The antenna 

C. L. Brown: The treatment of Gnenntionitiest in bonne deflection problems 

W. Z. Chien: The intrinsic theory of thin shells and plates . 

and A. Weinstein: On the vibrations of a clamped plate under tension 

L. H. Donnell: A chart for plotting relations between variables over their en- 
tire real range. 

H. L. Dryden: A review of the statiatical 

W. Feller: On A. C. Aitken’s method of interpolation 

K. O. Friedrichs and J. J. Stoker: Forced vibrations of einnin with 9 non- 
linear restoring force . 

R. E. Gaskell: On moment in 

A. Gelbart: (See L. Bers) 

G. H. Handelman: A variational principle for a state of combined plastic stress 

G. E. Hay and W. wets On vanes rigid frames loaded ataiaeies to 

M. Herzberger: A direct image error theme 

N. J. Hoff: A strain energy of derivation of the bedding loads 
of straight columns of thin-walled open sections 

Th. v. Karman: Tooling up mathematics for engineering. 

P. W. Ketchum: On the discontinuous flow around an airfoil wich hos 

W. M. Kincaid and V. Morkovin: An application of orthogonal moments to 
problems in statically indeterminate structures — 

Y. H. Kuo: On the force and moment acting on a body in deine how , ‘ 

J. Lehner and C. Mark: An cla of the method of the acceleration 
potential : 

C. C. Lin: On the motion of a ina fluid 

S. Lubkin and J. J. Stoker: Stability of columns and strings under oesladicaliy 

C. Mark: (See J. Lehner) 

V. Morkovin: On the deflection of anisotropic thin plates 

(See W. M. Kincaid) 

W. Prager: (See G. E. Hay) 

W. C. Randels: A new derivation of Munk’s formulae ‘ 

K. Riess: Electromagnetic waves in a bent pipe of rectangular cross soation. 

S. A. Schelkunoff: The impedance of a transverse wire in a rectangular wave 
guide 

- On the antenna problems 

H. J. Stewart: Periodic properties of the semi- permanent stmospherle peennae 


278 
281 | 
93 
168 
201 4 
349 j 
297 j 
61 
276 4 
7 q 
86 
97 
237 q 
351 
49 
69 
341 
2 
149 
334 
273 
250 
43 
215 
116 
328 
78 
354 
262 


J. J. Stoker: (See K. O. Friedrichs) 

— (See S. Lubkin) 

A. G. Strandhagen: Use of sine transform for non-simply supported beams 
J. L. Synge: On Herzberger’s direct method in geometrical optics . 

H. S. Tsien: Symmetrical Joukowsky airfoils in shear flow 

A. Weinstein: (See W. Z. Chien) 

Book reviews. 


346 
268 
130 


189 


| 
( 
‘ 
( 
{ 
| 
1 
=. 


281 


QUARTERLY OF APPLIED MATHEMATICS 


Vol. I JANUARY, 1944 No. 4 


THE TRANSFORMATION OF PARTIAL DIFFERENTIAL 
EQUATIONS* 


BY 


H. BATEMAN 
California Institute of Technology 


1. Introduction. In the early stages of the use of partial differential equations for 
the solution of problems of mechanics and physics the separation of variables and 
construction of simple solutions was the primary aim. The introduction of the idea 
of an exact differential by Fontaine and Euler led to the idea of associated differential 
equations such as those for the velocity potential and stream function in hydro- 
dynamics, the adjoint equations of Lagrange and Riemann, the contact transforma- 
tions of Legendre and Ampére, the transformations of Euler and Laplace for the 
solution of differential equations by definite integrals and other transformations too 
numerous to mention. Another aim which led to the study of transformations was 
that of reducing an equation to a canonical form. Laplace’s reduction of a linear par- 
tial differential equation of the second order to a form in which only one partial 
derivative of the second order occurs led to the study of transformations which pre- 
serve this form and of quantities which have a property of invariance. Conditions 
were then found that an equation may be reducible by means of a specified type of 
change of variables to some particular equations which had been fully studied. The 
conditions found by Campbell [1]+ (constancy of his two invariants) that Laplace’s 
canonical equation may be reducible to the equation of Euler and Poisson, may be 
cited as an example. 

A classification of transformations may be made by including in group A all 
transformations which arise from the condition or conditions that a linear differential 
form may be of a specified type (for example an exact differential). Transformations 
arising from the study of a number of linear differential forms may be included in this 
group. Transformations associated with the Calculus of Variations are also included 
because the equations of Euler and Lagrange are closely associated with the condi- 
tions for an exact differential. The extension of Legendre’s transformation found by 
Carathéodory [2] may be mentioned here. In this article attention will be devoted al- 
most entirely to transformations of group A. 

Transformations of group B include all those which arise from the conditions that 
a quadratic differential form may be of a specified type. The transformation of a 


* Received March 2, 1943. 
t Numbers in square brackets refer to the list of references at the end of the paper. 
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linear differential equation to a form in which the variables are separated is thus a B- 
transformation. The transformations of group B are not necessarily point transforma- 
tions, for instance, if Q(a, b,c) is a non-negative quadratic form in the real variables 
a, b, c transformations from (x, y, 2, t, u, v, w) to (X, Y, Z, T, U, V, W) may be con- 
sidered in which Q(dx—udt, dy ~~vdt, dz—wdt) goes over into Q(dX — Udt, dY — Vdt, 
dZ — Wat) the coefficients of Q in the first case being functions of x, y, z, t, u, v, w and 
in the second case functions of X, Y, Z, T, U, V, W. Since the equation Q =0 implies 
that dx =udt, dy =vdt, dz=wdt it also implies that dX = UdT, dY= VdT, dZ = WaT. 
In other words if u, v, w can be regarded as the component velocities of a recognizable 
moving particle of fluid then U, V, W can be regarded as component velocities of a 
recognizable particle of a corresponding fluid. Such a transformation is of interest 
because the density of each fluid can be defined in such a way that the equation of 
continuity is invariant under the transformation. 

Group C may be regarded as including all other transformations and some trans- 
formation of the other group which arise in the reduction of an equation to a canonical 
form. 

2. Associated equations of the types of Monge and Legendre. In his work on par- 
tial differential equations of the second order in two variables x, y which can be re- 
garded as independent, Monge [3] used z as dependent variable, p and gq as the first 
derivatives 2,, 2, respectively and 7, s, ¢ as the second derivatives 222, Zzy, Zyy. AS we 
shall have applications to fluid dynamics in mind, we shall deviate slightly from the 
notation of Monge and use u, v in place of p and g so that when 2 is the velocity po- 
tential u and v represent the component velocities as usual. This plan also allows us 
to use the symbol # to denote the pressure and g to denote the resultant velocity. 

The equations of steady motion of a compressible fluid under no body forces when 
the flow is irrotational and the fluid barotropic (density a function of pressure only) 
can, as we know, be derived from a variational principle 


sf p(u, v)dxdy = 0, Ss, 0 = Sy, (1) 


in which p is a specified function of g. For greater generality at the outset we shall 
suppose, however, that p is a specified function of u and v. Lagraiige’s partial differ- 
ential equation for this variational problem is then derivable from Haar's condi- 
tion [4] that p.dx—p,dy should be an exact differential. When the differentiations 
are made, the equation has the form 


Puut + 2puvs + Prout = 0. (2) 


When Legendre’s transformation is applied to this differential equation the new de- 
pendent variable is 

w=ux+ovy—2 (3) 
and since dz=udx-+udy, 

dw = xdu + ydv. (4) 


When wu and » can be regarded as independent this equation gives the relations 


Wy, y = (5) 
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and the equation for w is 
PuuWov — 2PurWuv + PovWuu = 0. (6) 


When, however, u and v are related so that they can be regarded as functions of a 
single variable r the equation (4) indicates that w is then also a function of r and we 
have the equations 


= xu(r) + yo(r) — 2, = + yo'(z), (7) 
which furnish a solution of (2) if 
+ + povd’*(7) = 0. (8) 
Let us now seek the conditions that 3 quantities 
R = R(u, 2), S = S(u, v), T = T(u, »), (9) 


may be the second derivatives Z,:, Zzy, Zz, of a single function Z(x, y). The required 
conditions R,y=S., S,=T, may be written in the form 


Rus + Rit =Surt+Sys, Sot = Tur + Ty. (10) 


We now seek the conditions that these two equations are both satisfied in virtue of 
equation (2). This will be the case when 


R, = R, = Wpuv + &, S,=- Whuuy S,=h-—- Wpur, (11) 
T, = Puus T, = Puv + S, = — Pury Pov. 
Equating the different expressions for Ry», Suv, Tus we obtain the equations 
hu = WuPuv — WePuuw hy = WuPpov — WePur, 
p p p p (12) 
hi = We Puu Pow hy We Puv We Pov. 
The elimination of h and h’ yields the two equations 
WuuPvy — 2WurPpuv + WovPuu = 9, 
p p p | (13) 
WuuPov 2WuvPur + WovPuu 0 


which show that w and w’ are solutions of equation (6). The case of chief hydrody- 
namical interest is that in which the second derivatives of w and w’ are all zero. We 
shall, however, look first at a possible alternative case. 

Equating the different expressions for S, and S, we obtain the equations 


f [ws d(p.) — d(p.)], 
(14) 
Pov = +f [ wud(p.) (pu) |. 


Hence 
Wy Pov + W = 2Wopuy + — 
Puv — We Puu + W = Wupuu + (15) 
Wu Pov + W Pure = Wepuu + 
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These equations lead to the relation 
wD, = wD,, where D = puupor — pro. (16) 
This equation is satisfied identically when D is constant but it may also be satisfied 


if w=E,, w’=E, where D and E are related. An important case of this second type 
occurs when D is a function of g only and w=v, w’ =u. In this case 


h=—-py S=—upy=— vp, (17) 
and # is a function of g only. If p. = —up, p» = —vp, where p is the density of the fluid 
we have 
R= — p — pv’, S = pur, T = — p — pu’, (18) 
and so Z is a kind of stress function which satisfies the equation 
RT — S? = p(p + pu? + pv*) = p(p + pq?) = F(R + T) (19) 
since R+T = —2p—pq?. The partial differential equation for Z is thus of Legendre’s 


type [5] 
5c(R, S, T) = 0. (20) 


In the special case in which — F(R+T7) = K?—}(R+T)? the equation reduces to 
one which occurs in Saint Venant’s theory of plastic bodies. This equation has been 
discussed by Hencky [6], Prandtl [7] and Carathéodory [8]. Oseen [9] uses the meth- 
od of Legendre in which the equation is first solved for R, differentiated with respect 


to y and so reduced to an equation 
(K? — V2)"(V.2 — Vyy) + 2V.V2y = 0 (21) 


in which 
T = Vy. 


It should be remarked that if p,=a#=2,, py= —2,, 


pb + pu? + po? = f, 1/p = —§, (22) 

we may write 
S=pa, T= — pa’, 
and the equations Ry=S,, S,=T, lead to the partial differential equation for the 


stream-function 2. 
In the theory of plane waves of finite amplitude equations of Legendre’s type oc- 


cur in at least two ways one of which is discussed by J. R. Wilton [10]. In the other 
way use is made of the equations 

R= op, S = — pu, T = pt pu’, (24) 
where now y denotes the time and x a co-ordinate in the direction in which the waves 
are travelling. The quantities u, » are again the derivatives of.a velocity potential 


2, p is the density, p the pressure and wu the velocity of the fluid. The quantities R, S, T 
are the second derivatives of a stress-function Z. The additional equations from which 


the relation between R, S and T may be derived, are 


v+3u? = f'(p), p= fle) — (25) 


- 
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The desired relation is thus 
T — RS? = f(R) — Rf’(R). (26) 
This equation, like that considered by Wilton, may be solved by the method of 
Legendre in which the equation is differentiated with respect to one of the independ- 
ent variables (in this case x) so as to reduce it to an equation of the Monge-Ampére 


type. The transformation to the new equation can be regarded as a special Backlund 
transformation [11] as Oseen [9] observes. If U=Z., the new equation is 


Uyy — 2(Uy/U2)U sy + (U2/U2) Use = — Rf" (R)U ce = CU (27) 


or 
HU,:+ 2KU.z+ LU, = 0, 


where H =(U;/U2)—c?, K= —(U,/U,), L=1. The invariant G is 
(28) 


and the condition G0 is satisfied so long as c?¥0. 
In the present case 


R, = — pu/c?, R, = — p/c*, Su = p(u?/c? — 1)S, = pu/c?, 


= pu(l — u*/c?), T, = — p(1 + u?/c’). 
The two equations (10) are both equivalent to 
(u? — c?)r + 2us +t =0, (30) 
and it is readily seen that w= —1, w’ = —u, g=k?—hl =c? 0. 
It should be noticed that if we solve equations (9) for u and v in the form 
u = F(R, S), = G(S, T), (31) 
the equation u,=v, is satisfied on account of Ry=S., Sy=T; if 
Fr =Gs, Fs =Gr. (32) 


These two equations then are consequences of the single equation #(R, S, T) =0. The 
expression of such an equation in the two forms (32) may be regarded as a problem of 


some interest. 
In the case when D is a constant and » is a function of g only 
D = = (q/c)?], (33) 


and the flow is either entirely subsonic (D>0) or entirely supersonic (D <0). In 
many cases in which is a function of g only, D can have either sign and so the flow 
is partly subsonic and partly supersonic. It is then of some interest to seek the condi- 
tion satisfied by the function 3¢(R, S, T) when D>0. For this purpose we write the 


equation in the form 
0= K(R,S, T) = [(R — T)? + 482]? — J(R+ T), (34) 


where J is a function which is such that 


2ef'(e) = — J[— 2f(p)]. (35) 
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Here »=f(p) —pf’(p) is the relation between the pressure p and density p. Now we 
find by differentiation that 


+ Rr —2= —2—2'(R+T) = — 4 — 2f"(0)/f'(o) = 4(c2/q? — 1). (36) 
Hence Ke+Hr—2>0 when c?>g? and <0 when c? <q’. In the case of the 
plastic equation J is a constant and so 

Re + Hr = 0. (37) 
The corresponding flow is characterized by the relation g?= 2c? and is consequently 


supersonic. 
A simple case in which D is constant is obtained by writing 


p = au? + 2cuv + bo’, (38) 
where a, 6 and ¢ are constants, The functions w, w’ both satisfy 
bwWun + AW, — = O, (39) 
and we may write w=)dV,, w’=aV,, where V is a solution of this equation. If a func- 
tion W is defined by the equations 
W, = cV, — aV,, W, = bV. — (40) 
we may write =2)W, and it is readily found that we can write 
R= 20(W + cV), S = — 2abV, T = 2a(W + cV), (41) 
where V and W are connected by the foregoing equations. In this case the relation 
between R, S and T is simply 
ax = aR — bT =0. (42) 
The quantity e+X7—2 is now simply —(a+b)/a, a constant. There is no change in 
sign of the expression. It will be noticed that the equations aR —bT =0 and 
ar + 2cs + Dt = 0 (43) 
satisfy the condition of apolarity 
Ab + Ba — 2Cc = 0, (44) 
when the first equation is written in the form AR+2CS+BT=0. 
3. The transformation of the Monge-Ampére equation. If for the equation 
hr + 2ks+lti+m-+ n(rt — s*) = 0, (45) 
g= k? —hit+tmn (46) 


the expression 


is not zero and so the two systems in the methods of Monge and Boole are distinct, 
the equation is transformed by a contact transformation 


X = X(x, y, 2), Y = Y(z, 9, #, 2), Z =Z(x, ¥, (47) 
U = U(x, y, 2, u,v), V = V(x, y, 2, u, 2), dZ =UdX —VdY = o(dz — udx —vdy) 
into an equation 

HR+2KS+LT7+M + N(RT — S*) =0 (48) 


} 
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for which the quantity G=K*—HL+MN=0. 

In a paper published in 1904 Sophus Lie [12] remarked that it would be desirable 
to have a direct proof of this theorem and Kiirschék [13] gave one based upon a 
representation of the equation in the form of a Jacobian 


d(a, b)/d(x, y) (49) 


where a and 6 are functions of x, y, 2, u, v and d/dx=0/0x+4u(0/0z) +r(0/du) 
+5s(0/dv), d/dy=0/dy+v(0/dz) +5(0/du)+t(0/dv). When this representation is not 
used the proof is algebraically more difficult but the analysis is worth giving on ac- 
count of the numerous relations to which it leads. Reference for this type of proof may 
be made to a paper by R. Garnier, Sur la transformation des dérivées secondes dans les 
transformations de contact et les transformations ponctuelles, Bull. des Sci. Math. (2), 
64, 13-32 (1940). 

We shall suppose that dz = udx+vdy and that consequently dZ = UdX + VdY. To 
make dU=RdX+SdY, dV=SdX+TdY consequences of du=rdx+sdy, dv=sdx 
+tdy we shall require that 


dU — RdX — SdY = (Uy — RXu — SY«)(du — rdx — sdy) 
+ (U, — RX, — SY.)(dv — sdx — tdy), 


dV — SdX — TdY = — SXu — TY.)(rdu — dx — sdy) 
+ (V, — SX, — TY,)(dv — sdx — tdy). 
With the notation 
Z,:22,+ a2, Z. = Z, + etc., 
(uu) = Uy — RX, — SYu, (uv) = U, — RX, — SY», (51) 
(ou) = Va — — TY x, (w) = V, — SX, — TY,, 
the equations to be satisfied are 
r(uu) + s(uv) + — RX, — SY, =0, r(uv) + s(vv) + Vi — — TY: = 0, 
s(uu) + t(uv) + U2 — RX2— TY, =0, s(vu) + t(vv) + V2 — SX2 — TY 2= 0. 
Hence 
2), 
tA= UD u—V2U ua tR(XwW2— u) 
+T(U.Y1—UiY .) 


sA=V2U,—UWV V2) +S(U2X VV 
A(rt—s?) = UV 2— + T(U2Y1— 2) 
1) 
+5S(V u— 0) +(RT—S*)(XuV Xe¥ u)- } 
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i 
| 
‘ 5 
| 


288 H. BATEMAN [Vol. I, No. 4 


The two expressions for s are equivalent on account of the relations 
[uv] = [xv] = [yu] = [xy]=0, [xU]=[YV]=c (54) 
which, in addition to the relations [YZ] =[ZX]=0, [UZ] =cU[VZ] =cV are satis- 


fied because the transformation is a contact transformation. In these relations [AB] 
is the Poisson bracket 


[AB] = — + AvB2 — A2B,. (55) 


The relations are derived by Lie [12] by a clever device. In the book of Cerf [14] the 
relations are derived from the equation 


o[AB] = [adj (56) 


where A, B are the expressions for a(x, y, 2, p, g), b(x, y, 2, p, g) in the new co-ordinates 
[X, Y, Z, P,Q], while F. Engel [15] obtained them with the aid of the bilinear co- 
variant by a development of a method used by G. Darboux. 

It is readily seen that the equation hr+2ks+li+m-+n(rt—s*) =0 becomes HR+2KS 


+LT7+M+N(RT—S)=0, where 
H = hP, + k(Po + + + mP, + nP., 
2K = h(Rp+Cp) + k(Co+ Ro +C, + Ry) + + Ca) 
+ m(R, + C,) + n(R. + C.), 
L = + + Aq) +144 + mA, + nA,, 
M = hQ, + k(Qs + Qa) + Wa + + nQ., 


(57) 


N = hB, + + + 1B. + mB, + nB., 
where 
Ag = — U2Yu, Ay = — A, = U2¥,; — 
A, = UiY, — UY, A, = U2:Y, — U.¥2, A, = — 


— X2V1, 


B, = X.Y: — B, = X2¥», 
Poe PeXVi-XVe, XVi— 

24, P, = XV, — XV2, 


U2, — Q, = UiV, — Ui, Q. = — U2Vi, 


> 


Q> = ULV; = UV 2 U2)», 0, = UV «, 
YuVe2 — YoV iu, Ry = — ViVi, R. = Y2V1 — 
Y2V, — YWVu — o. 


R, = — YVi, R, 


These equations give the relation 


ta 2 
= 

far: 
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K*—HL+MN=J(k?—hi+mn), if ) 

(Pot Py) (AstAg) =J, 
3(R-+C,)(Re+C.) + (Q,Be+O-Br) 
PcAa=9, 
(Co+Rot+C,+R,) +0,(B.+B,) +B,(Qo+Q,) 

= P,(Ast+A,) +A Py), 
+B.(Q.+Q,) 

= 
+ C,)(Co+Rs+C,+ Ry) +O, (Bo+ By) +B, 

= 
+O. (Bo+Bz) +B. 

= 
PpA,— P,A,=0, 
3(Rp+C,)(R-+C.) +0,B.+0-By— PpAc— PA y=9, 
PeAc—P-Aa=0. 


These relations may be established by using a parametric representation of the 
quantities satisfying Lie’s conditions for a contact transformation, we therefore write 


X, = aye + are’, X2 = bie + dee’, Xu = + = die + dee’, 
aif + def’, U2= bif + bef’, U.= af + cof’, U,= dif + df’, 
— VY. = asp + ap’, = bsp + bap’, Yi=csp tsp’, Vo = dsp + dap’, 
— Vi = ag + aq’, + bq’, Vi=esg V2=dsq + dag’, 


where the quantities 


(58) 


(59) 


G2 
b, be bs by 
C1 Co C3 
d, dz dy 

form an orthogonal matrix and the quantities e, e’, f, f’, p, p’, g, g’ are such that 
— Coa, + — (ef’ + e’f) = (asa4 — + — buds) (pq’ — p’g) =. (60) 

It is then found that 
— H = eq(1, 3) + e’g(2, 3) + eg’(1, 4) + e’g'(2, 4), 
— L= fp(, 3) + f’p(2, 3) + fp'(1, 4) + 4), 
2K = (ef’ — e’f)(1, 2) + (pq — p’9)(3, 4), (61) 
M = fq(1, 3) + 3) + fa'(1, 4) + f'9'(2, 4), 

N = ep(i, 3) + e’p(2, 3) + ep’(1, 4) + e’p’(2, 4), | 


| 
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where 
(1,3) = h(dics + ayb3) — + C1d3) + R(did3 + bibs — — C1€3) 
+ n(aids — bycs3) + m(dya3 — cbs), 
(2, 3) = + obs) — + Code) + R(deds + bobs — a2a3 — 
+ — bees) + m(d2a3 — C2bs), 
(1,4) = h(dicg + ayb4) — + crds) + R(dids — byb4 — — 
+ n(aidy — bica) + m(dia, — C154), 


(2,4) = + Gods) — + Coda) + + baby — 204 — C2C4) 
+ n(a2d, — boca) + m(d2a4 — C2b4), 
cs, 2) = h(d,a2 d2a1) + 1(bi¢2 be¢1) + — + — d2b,) 
+ n(ayb2 + a2b1) + m(cid2 — 
(3, 4) = A(bscg — bacs) + — agdy) + R(csaq — + — 
+ n(dscq — dycs) + m(bsa4 — yas). 
It is readily seen that 
MN — HL = [(1, 3)(2, 4) — (2, 3)(1, 4) (ef — ef’)(0'¢ — 9’) (63) 
and that, on account of the properties of an orthogonal matrix 
— ef’)? = (p’q — pq’)? (64) 


The expression for K also simplifies considerably and the proof may be readily com- 
pleted. The quantity J as in Kiirschék’s analysis, is equal to o* and so is not zero. 

Contact transformations are not the only ones in which the condition g+0 is 
invariant. In the theory of the steady two-dimensional motion of an inviscid elastic 
fluid the equations satisfied by the velocity potential z and stream-function Z are 
respectively 


+ 2Puvs + Pout = 0, + 2Purs + Pot = 0. (65) 
In this case dz=udx+vdy, dy =p,dx —p,dy and so 
= po i= — pu; = — PuuPov = C(u? + v? — (66) 


Thus g=0 either when c=0 or when g=c. The supersonic region is characterized by 
the condition g>0 and the subsonic region by the condition g <0. The curve for which 
c=0 is a boundary for the flow just as in the case of the Prandtl-Meyer flow round a 
corner. The transformation under consideration is a special Backlund transformation 
and is included in the group of Backlund transformations 


= X(x, y, u,v), J= V(x, y,u,v), &#= U(x, y, u,v), = V(x, y, u, 2), (67) 


for which the Jacobian 0(X, Y, U, V)/0(x, y, u, v) is not zero. These transformations 
have been studied carefully by Goursat [16]. The requirement that 


should be exact leads to an equation of the Monge-Ampére type in which z does not 
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occur explicitly. It is shown, however, that the general Monge-Ampére equation of 
this type cannot be obtained in this way and a similar result has been found by J. 
Clairin [17] in his studies of more general Backlund transformations. Clairin has 
studied in particular transformations of type 


x’ = fi(x, y, 0,2’), = fo(x, y, 2, 2, (68) 
u’ = fs(x, 2, U, 2’), x, y, 2, 2’). 


Some of Clairin’s work is summarized in the book of Forsyth [18] and illustrated by 


means of examples. 
Another transformation of type (67) which preserves the condition g¥0 is ob- 


tained by writing 
dU = Rdx + Sdy, dV = Sdx + Tdy, 
where R, S, T are the functions of u and v used in section 2. Making use of the equa- 
tions 
(upu + vp.)dx = pudz + (upu + vp.)dy = prdz — udz, (69) 
we find that 


(Tu — Sv)dU + (Ro — Su)dV = (RT — S*)dz, | (70) 


(Spu + Tp.)dU — (Rp. + Sp.)dV = (RT — S*)di. 


Hence, if 
=U, y'=V, =02/8U, v =02/dV, = = 92/8V, (71) 


we have the relations 
u’ = (Tu — Sv)/(RT — S*), vy = (Rv — Su)/(RT — S?), 
a’ = + Tp.)/(RT — S?), (Rpu + Sp.)/(RT — S?) 


which, in conjunction with the preceding relations define two Backlund transforma- 
tions. The transformations considered in my paper on the lift and drag functions are 
of this type [19] and are not generally contact transformations as is apparently im- 
plied by a statement relating to the correspondence of the supersonic regions in the 
two associated types of flow. 

In the case in which p is a function of g only the relations between wu’, v’, u, v are 


(72) 


u=—u/p, v=—r/p, 
and, if p’= —1/p, p’+p’g” = —1/(p+pq*) we have 
+ 
dp’/dq’ = — qp/(p + = — pg’, 
1 — q’2/c’? = (1 — q?/c*)[p/(p + 9%) where = dp'/dp’. 


This transformation may be compared with that obtained by means of Haar's 
adjoint variation problems [20]. In this case 


u* = pu/(p — — Vpy), v= UP»), 1/(p — — Vp»), 
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and, when p depends only on gq, 
u* = — pu/(p + pq’), v* = — pv/(p + pq’), 
p* = + 1/(p + 9°), q* = q/(b + 9). 
Defining p* by the equation dp* = —g*dq*, we have 
= (p + 0q*)/pp, = p* + = 1/f, 
ct? = p2pg?(g? — c*)/(p + pq*)?[c*(p + q*)? + p2(q? — c*)], 
— g*? = — + + p*(q? — c*)]. 
Hence c*?=g** when c?=0 and c*?=0 when g?=c?. It should be noticed, however, 


that 
— g?) = — g*)p*ptg*/(p + + + — 


This may be compared with Haar’s general relation* 


Transformations more general than those of Bicklund have been considered by 
Gau [21] but so far no hydrodynamical applications have been found for these so far 
as I know. Mention should be made, however, of the equiareal transformations from 
the Eulerian to Lagrangian co-ordinates in the two-dimensional flow of an incompres- 
sible fluid. These transformations have been much used in mapping but the hydro- 
dynamical applications are beset with formidable difficulties. 

No mention has been made of the use of transformations in the theory of surfaces, 
congruences, etc. This is a subject which is well treated in the books of Darboux [22], 
Forsyth [18], Goursat [16], Bianchi [23] and Eisenhart [24]. 

4. Transformation of the linear equation. In the special case in which A, k and | 
are functions of x and y only, n=0 and m is a linear homogeneous function of u, v 
and z with coefficients depending only on x and y, the Monge-Ampére equation reduces 
to a linear equation. The behavior of this equation in transformations of type 


X = X(x, y), Y = Y(zx, y), Z = 2F (x, (73) 


has been studied by Darboux [22], Cotton [25], Rivereau [26], J. E. Campbell after 
the case F=1 had been discussed by Laplace [27], Chini [28] and others [16]. 
Campbell uses the equation in a form in which g = 1, a form to which the general equa- 
tion can be reduced by multiplying it by a suitable factor. He then shows that there 
are two invariants J, J and an absolute invariant J/J where if suffixes denote partial 


derivatives 


ha, + k(ay +b.) + lby + (hz + kya + (kz +1,)b 

+ ha® + 2kab + 1b? — m,, J =a, — bz, (74) 
2a = + ky — my) — +1, — m,), 
2b = +1, — m,) — R(hz + hy — mz). 

* This is a consequence of the relations 


v) = pp (PuuPoo Pur); a(u*, v*) p PU Po Pu v ). 


In the correspondence between the two hodograph planes complications arise on account of the relation 
Of PuwPus— to these Jacobians. 
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Laplace's invariants are }(I—J), }(1+J). These are used with a different notation, 
in Darboux’s Théorie des Surfaces, t.2. Campbell shows that in the case of the equa- 
tion of Euler and Poisson the invariants J and J are constant. This may be compared 
with Cotton's result. The harmonic equations belong to the group characterized by 
the relation J=0. The equations considered by Burgatti [29] are such that J=0. 

In mathematical physics the simple solutions of linear equations play an impor- 
tant part and the primary problem is that of separability. Even in the case of the 
equation with two independent variables there are some unsolved problems. A good 
idea of the progress which has been made may be derived from Darboux’s book [22]. 
The method of Laplace provides an important way of reducing equations by a cas- 
cade process which is particularly useful in the treatment of equations arising in the 
theory of plane waves of finite amplitude. Reference may be made to a paper of 
Love and Pidduck [30], an article by Platrier [31], some papers by Bechert [32] and 
to two papers by Oseen [9] in which the transformation and reduction is given for 
equations occurring in the theory of earth pressure and in the theory of plasticity. 

In the theory of the steady motion of an inviscid compressible fluid the equations 
in the hodograph plane are linear. These equations are 


PovWuu — 2PuvWur + PuwWov = 0, Povtaa — 2pwta + pula = 0, (75) 
= UW, + VW, — W= — z= twa + 105 — W = — W. 
When # is a function of g only the equations #=p,, v= —p, take the form 
i = — pt = — Gsinr, = pu = cost, = pq, (76) 
and the equations become 
Wrr + = + q(qwz), = 0 (g a function of q). (77) 
These are consequences of simple relations between w and # 
Ww, + gw; = 0, W, — qw; = 0. (78) 
The corresponding relations between z and are found to be 
= = — Gar (79) 
and so the equations for z and 2 are 
(q°/4) + = 0, 
(g°/9) + = 0. J 


These are equivalent to the equations obtained by Molenbroek [33] and Tschaplygin 
[34] for the case in which the relation between p and p is of the polytopic or adiabatic 
type. The symmetrical forms of the equations are easy to remember. 

It is sometimes useful to introduce other quantities which satisfy linear relations. 
Thus we may obtain the desired relations between 2, 2, w, by writing 


(80) 


w= = W=e, = — 


—(q/Gerr, 2= (G/qerr — ery 


2 


where 
Crr + G(qea)q 0, err + q(Geq) = 0. 
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The literature dealing with the transformation of linear equations in several vari- 
ables is very extensive and only a brief summary can be attempted here. Beltrami's 
work on differential parameters [35] was extended by Ricci and Levi-Civita [36], 
Cotton [25], Levi-Civita [37] and many other writers. The development of general 
relativity, electrodynamics and the theory of elasticity has made this work more or 
less known. The work of Lamé on simple solutions of the potential equation [38] was 
much developed by later writers and a good summary of results up to 1893 is given in 
the book of Bécher [39]. The use of a variational principle for obtaining the trans- 
formation of the equation was recommended by Larmor [40], Volterra and others 
[41]. Since the advent of the new quantum theory the interest in separable equations 
and separable systems has much increased. Mention may be made of the work of 
Staeckel [42], Eisenhart [43] and Robertson [44]. 

In addition to the simple solutions of partial differential equations there are solu- 
tions having the form of products in which one or more or the factors satisfies a par- 
tial differential equation instead of an ordinary differential equation. Comparatively 
little work has been done on this problem. In the case of Laplace’s equation Vzz+ Vyy 
+ V,.=0, the aim is to find a solution of form [45] 


V = ZF(X,Y), (generalized binary potential) 
where F satisfies a partial differential equation of the second order in the variables 
X and Y. The problem seems to depend on the formation of a relation of type 

(p? + g? + 1?) (dx? + dy? + — (pdx + qdy + rdz)? = adX*? + 2hdXdY + bdY? 
in which a, } and h are functions of X and Y only. There is a similar relation for the 


corresponding problem in any number of variables. 
5. The transformation of Legendre’s equation. Legendre’s equation 


S, T) = 0 
is unaltered in form by a Legendre contact transformation 
X' = U, U' = X, Z’= UX + VY --Z, 
which makes 
R’ = T/(RT — S®), S’ = — S/(RT — S?), T’ = R/(RT — S?). 
In particular, the equation 
R'T’ — S" = F(R’ + T’) 
R+T 
RT — 3) 
an equation of the same general type. Again, if a, b and h are constants the contact 
transformation 
Z’ X’=X, 
U' = aX + + U, V’ = 


makes R’=a+R, S’=h+S, T’=b+T and so transforms an equation of Legendre’s 
type into another equation of the same type. Equations of the preceding type 
usually go into Legendre equations of a slightly different type. 


becomes 


1/(RT — S?) - F( 


| 
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Other transformations may be found by first transforming the equation to the 
Monge-Ampére form by Legendre’s device. If, for instance, the equation is 


T = F(R, S) 
and we differentiate with respect to x using then the new notation z= U, R= U,=u, 
S= U,=v0, T,=S,=t, Rz=r, S:=s, the new equation is 
t= F,,(u, + F,(u, v)s. 
Comparing this with the equation +2puvs + hoot =0 we find that 
Fy = — pus/Pov, Fo = — 2pus/Pov. 
Eliminating F we find that the equation for z is not a general equation of the type 
considered in §2 because the function p(u, v) satisfies the condition 
= 0 where D= — 
An equation for which D is constant satisfies this condition and the equation of 
type 
a(R, S, T) = 0 
associated with it by the analysis of section 2 may be regarded as a transform of the 
original equation 3(R, S, T) =0. In the case when D=1, we may write 
fuu = e* sec Pov = sec Pu» = tan 
where a and # are functions of u and v which must be chosen so that 
sec? bb, = e* sec b(a,) + e* sec b tan b(b,), 
sec? bb, = — e~* sec b(a,) + e~* sec b tan b(b,). 
Also, since F,, = —e**, F, = —2e* sin v, we must have the additional equation 
2e* cos b(bu) — 2e7*(a,) + 2e* sin v(a,) = 0 
which is seen, however, to be a consequence of the other two. Elimination of the 
derivatives of a gives the equation 
Cb yy + — 2 sin =O OFr + Povbuu — = 
and it is readily seen that a satisfies the same equation. The equation D=1 is given 
as an example in Forsyth’s book, p. 220, Ex. 11. 
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THE INTRINSIC THEORY OF THIN SHELLS AND PLATES* 
PART I.—GENERAL THEORY 


BY 


WEI-ZANG CHIEN 
Department of Applied Mathematics, University of Toronto 


1. Introduction. The development of the theory of thin shells and plates by many 
authors [1, 2, 3, 4] can be summarized under the following three main heads: 

(1) All theories are based upon certain simplifying unproved assumptions. For 
example: (a) the thickness remains unchanged during the deformation, (b) the middle 
surface in the unstrained state deforms into the middle surface in the strained state, 
(c) the normals of the unstrained middle surface deform into the normals of the 
strained middle surface. 

(2) All theories involve the use of displacement to describe the state of deforma- 
tion. This plan works well in the theory of small deflection, but presents considerable 
difficulty in the case of large deflection. 

(3) The various approximations used in the theory of thin shells and plates are 
confusing. If one attempts to give a complete picture of the theory, one must be able 
to introduce a systematic method of approximation, which not only clears away the 
confusion of various approximations, but also leads to a complete classification of all 
thin shell and plate problems. 

The purpose of this paper is to give a systematic treatment of the general problem 
of the thin shell, which includes the problem of the thin plate as a special case. The 
work is based on the usual equations of elasticity for a finite body, supposed to be 
homogeneous and isotropic. The final equations of Part I are the three equations of 
equilibrium, (6.34), (6.35), and the three equations of compatibility, (6.43), (6.44), 
for the six unknowns, pas, Gas, Which represent extension and change of curvature of 
the middle surface. When these quantities are found, the strain and stress throughout 
the shell or plate can be calculated. The displacement does not appear explicitly in 
the argument. Since we deal rather with stress, strain and curvature (all tensors), the 
tensor notation proves much more convenient than any other. 

In Part II and III, we shall discuss the various approximate forms of the equations 
arising from consideration of the thinness of the shell or plate and the smallness (or 
vanishing) of its curvature. The strain is, of course, always supposed to be small. We 
obtain a complete classification of all shell and plate problems. There are found to be 
twelve types of plate problems, and thirty-five types of shell problems. To each type, 
there corresponds a set of six equations which are simplifications of the equations 
(6.34), (6.35), (6.43), (6.44) of Part I, with certain terms dropped on account of 
smallness and the uncalculated residual terms omitted. The equations obtained in- 
clude all the familiar equations in the field of small deflection and the few equations 


* Received May 13, 1943. This paper is part of a thesis submitted in conformity with the require- 
ments for the degree of doctor of philosophy in the University of Toronto, 1942. 
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already known for large deflection. The new results for finite deflection may prove 
particularly interesting. 

The author is much indebted to Prof. J. L. Synge who not merely has directed this 
work but has actively participated in it. 

2. Reduction of force system to reference surface; macroscopic equations of 
equilibrium. We shall start in this section by reviewing some main results in a previ- 
ous paper [5]. 

All theories of thin shells and plates involve the use of a reference surface. Usually 
the middle surface is taken without explicit distinction between the middle surface 
in the unstrained state and the middle surface in the strained state. In the present 
methodical treatment we shall use a general reference surface in the material in sec- 
tions 2-5. In section 6 and later parts we shall use the surface in the strained state 
formed by the particles of the middle surface in the unstrained state. To the order 
of approximation used there, this is indistinguishable from the middle surface in the 
strained state. (This is generally assumed to be the case; cf. E. Reissner [4].) 

The following coordinate system will be used: x° at any point A inside the material 
of the shell or plate is the perpendicular distance of A from the reference surface Sp, 
and x* at A are the values of any Gaussian coordinates on Sp corresponding to the foot 
of the perpendicular dropped from A. (Throughout the paper, Latin indices have the 
range 0, 1, 2 and Greek indices the range 1, 2; summation over either of these ranges 
will be signified by the repetition of an index.) This may be called a normal space co- 
ordinate system with respect to Sp. 

Let us denote the line element in space by ds?=g;;dx‘dx’, where gi; is the funda- 
mental tensor. Furthermore let gjo)i; be the values of gi; at So, then we have in usual 
tensor notation 


Z00 = = = = 1, Soa = = = = QO, (2.1) 
= Gag, = 2%, det. = gto) = det. (dag) = a, (2.2) 


where d,3dx%dx* is the metric on Spo. 

We shall now consider forces in the shell or plate. Let C be a curve on Sp and Ag 
a point on C. Let mj) be the unit vector in Sp normal to C at Ag, indicating the posi- 
tive side; let Af) be an arbitrary unit vector in So at Ao. We consider the system of 
forces acting across an element of area standing on the element dso of C at Ao, and 
terminated by the surfaces of the shell or plate. We replace the forces acting on the 
element by a statically equipollent system acting at Ao. This leads to the following 
invariants, which in fact define the macroscopic tensors, TJ*° (shearing stress tensor), 
T*8 (membrane stress tensor), L** (bending moment tensor): 


T n1o}2d59 = shearing force normal to So across the element dso, 
= component in the direction of of membrane stress across 
the element dso, (2.3) 


component in the direction of Afo, of the bending moment 


across the element dsp. 


Similarly, let dS, be an element of area of So at Ao. Consider the external forces 
acting on the volume element consisting the normals to S» standing on the element 
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dS, and terminated by the surfaces of the shell or plate. We replace them by a stati- 
cally equipollent system acting at Ao. This leads to the following invariants, which 
give the definitions of the external force and external moment tensors, Fi, M*: 


F°dS) = normal component of the external force on dSo, 
FA .o}4@59 = component in the direction of Af, of the external force on dSo, (2.4) 
MA (o}4dS9 = component in the direction of A{o, of the external moment on dS». 


Then from purely statical considerations, we obtain the following six equations of 
statical equilibrium ([5], p. 110): 


(a) Fo = 0, 
(b) + + F* = 0, (2.5) 
(c) + + Me = 0, 
(d) ntojasl — = 0. 
The symbols have the following meanings: 
Qa 
= 7”. { \ 
Ta) 
a B 
= + { { T@*, (2.6) 
TA) a TA} a 
bas = 2%=0, 
Nlojap = €11 = €29 = 0, = — = 1. 


The Christoffel symbols are calculated for aas. The quantities (1/2)bas are the coeffi- 
cients of the second fundamental form of So; they vanish if Sp is a plane. The radius 
of curvature R in the direction of a unit vector yj) (counted positive when Sp is con- 
vex in the sense of x° increasing) is given by 


2/R = bas (2. 7) 
By eliminating T** from (2.5) we get a set of three equations, 


+ + F* + = 0, (2.8a) 
— T™ + + = 0. (2.8b) 


These equations, rather than the equations (2.5), are fundamental in the later theory. 
It should be noted that in the case of repeated covariant differentiation with respect 
to das, the order of the operations cannot be changed unless the total curvature K 
of the reference surface So (cf. Eq. (3.13)) is equal to zero. 

The above equations are valid for shells and plates of finite thickness. When we 
come to deal with approximations based on the smallness of certain quantities, we 
must of course consider only the magnitudes of dimensionless quantities. It is 
best therefore to work with dimensionless quantities throughout. Let us introduce a 
standard length L, a lateral dimension of the shell or plate (e.g. the diameter in the 
case of a circular plate). By ds we shall understand the distance between two adjacent 
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points, divided by L; this dimensionless ds may be called the reduced distance. Simi- 
larly all coordinates x‘ are supposed to be in reduced or dimensionless form. Conse- 
quently the fundamental tensor g;; is dimensionless, and all tensor operations (such 
as raising or lowering suffixes or covariant differentiation) are dimensionless opera- 
tions. We also reduce stress to dimensionless form by dividing by Young’s modulus £, 
and body forces by dividing by E/L. 

All the relations written above hold equally well for reduced or dimensionless 
quantities. We shall carry through the work with these quantities; if we wish to trans- 
late conclusions into ordinary dimensional form, we have simply to multiply by that 
combination of the form L”E£” which restores the required dimensionality. Young's 
modulus will not appear explicitly in our work, since it is eliminated from the stress- 
strain relations by the process of reduction. 

We also note that a thin shell or plate is defined as one whose thickness is small in 
comparison with a lateral dimension LZ. Customarily, a thin shell is defined as one 
whose thickness is small compared with the radius of curvature; this is unsatisfactory 
in the limiting case of a plate, and also in the case of a shell whose thickness is small 
in comparison with the radius of curvature but of the same order as the lateral di- 
mension L. 

3. Representation of 7°°, T°, L**, Fi, M« as power series in the thickness. Let C 
be a curve on the reference surface So, and 2 the surface formed by erecting normals 
to Sy along C. Let dso be an element of C, and d= 
the strip formed by the normals on dso (Fig. 1). 
Let ds be the length of the arc of intersection of d= 
and the surface x°=constant, passing through any 
point A in dd. Let nu, be the unit vector normal to 
d> at A, and m,o)4 the unit vector normal to d= at 
the reference surface So. 

Here we note that there are two distinct classes 
of quantities: (1) those defined only on the refer- 
ence surface, such as Gag, bas, T*°, 
Fi, Ej, (normal derivatives of the stress tensor 

Fic. 1. on So); (2) those defined at any point in the ma- 

terial of the shell, such as g;;, E;; (stress tensor), 

X; (parallel vector field), 2. For all space tensors, the change of indices will be effected 

by applying g;; or g‘7; and for all surface tensors, the change of indices will be effected 
by applying gjo);; and gi, [defined as in (2.1), (2.2)]. 

By the definition (2.3), the shearing stress tensor and the membrane stress tensor 
are calculated by the invariant formula: 


T2'n ad [0] i459 = f (3.1) 
h 


where X; is any parallel field of unit vectors, X10}; the same vector field at the reference 
surface So, and E*‘ part of the stress tensor E‘/. The symbol h under the sign of in- 
tegration indicates here (and throughout the paper) that the integral runs from 
x°=—h_) to x°=+h,,), where both and are positive functions of x* (for 
thin shells or plates, they are assumed to be small). Furthermore, the bending moment 
tensor is caculated from the invariant formula: 
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h 


where 


= g = det. (gas), 


ell = 22 = Q, eit — of! = I, 


(3.3) 


For the external force system, assume that X#), is 
the component of the body for per unit volume in the 
arbitrary direction of \,; at any point in the shell, and 
Zix) and Z{_) the components of the given 
loads per unit are applied to the upper and lower sur- 
faces respectively. Let us consider (Fig. 2) a portion of 
the shell or plate obtained by drawing normals to its 
reference surface over a surface element dSo. The por- 
tions of the boundary surfaces of the shell or plate cut 
out by these normals are do,,) and do,_). Let dS be the Fic. 2. 
corresponding element drawn at constant normal dis- 
tance from the reference surface Sp. Then the external force and moment com- 
ponents are calculated from the invariant formulae: 


(0) :dSo = f + Zi hy + (3.4) 
h 


M = f (npad*X8dS) + 
h 


where (4,8, %(—)as are the values of nag at the upper and lower surfaces of the shell 
or plate respectively. 
In order to carry out the integrations in (3.1)—(3.5), we must express all the quanti- 
ties in the integrands as functions of x°. These can be written as follows: 
Sas = Gag + bapx® + 3Cap(x°)?, (3.6) 
a{1 + 2Hx® + K(x®°)?}2, 


(3.7) 
d 
Na = = 1 + + K(x*)*}, (3.8) 
So 
Ao = Apoyo, = + (3.9) 


The first relation is well known; dag, (1/2)bas, (1/2)¢ag are the first, second and third 
fundamental tensors of the reference surface So respectively. These tensors are not 
independent, but connected by six geometrical conditions of flat space ([5], p. 112). 


Three of these read 
Cas = (3.10) 


and the other three are the well known equations of Codazzi and Gauss, 
baris = 0, (3.11) 
AR, apy = Dpsbay — Opybas. (3.12) 
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Here the single stroke indicates covariant differentiation with respect to x* and the 
tensor dag; Ryasy is the two dimensional curvature tensor formed from the tensor dag 
(sometimes called the Riemann Christoffel tensor of the surface So, see [6], p. 182, 
Eq. (50)). 

The relation (3.7) can be obtained by direct calculation from the definition of g. 
Here K is the total curvature and H the mean curvature of the reference surface So. 
If Ri, Re are the principal radii of curvature, then 


1/1 1 1 
H= K= (3. 13a) 


therefore, in tensor notation, 
SK = 47 = (3.13b) 


We also note the following relations, which are often used in later calculations, 


(a) = aa”, 
(b) = — 58), (3.14) 
(c) = 2a{ (4H? — K)a® — Hb}, 


The proof of the geometrical relation (3.8) is long, but not difficult. The relation (3.9) 
is obvious; since \,; is a parallel vector field, we have 


Agi = 0. (3.15) 


Here the double stroke indicates covariant differentiation with respect to the space 
coordinates x‘ and the tensor g;,;. Putting 7=0, =a, 7=0 in (3.15), the relation (3.9) 


follows at once. 
Besides all the relations (3.6)—(3.9), we also need the following geometrical results: 


dS 1/2 da 1 do 1 
’ = ’ = = 0 (2) (3 16) 
dS a dSy |N%x,|\a dS |NL,|\e 
the positive root being understood. Here N,4)i, Nii are unit normal vectors, drawn 
out from the upper and lower boundary surfaces respectively, and g,+), gi—-) are the 


values of g on the boundary surfaces. 
Substituting (3.8), (3.9) in (3.1), we have 


1/2 
f + 3b§x°) (£) aah 
h 


1/2 
+ more f aah = 0. (3.17) 
h af 


But mo}2 and Ajo}; are arbitrary, and consequently we have for the shearing stress 


tensor 
h a 


and for the membrane stress tensor 
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1/2 
Tos = f + 468 2°) (3.19) 
A a 
Similarly, substituting (3.6)—(3.9) in (3.2), we obtain the bending moment tensor 
g 1/2 
= an f + 29 +) (3.20) 
h a 


Furthermore, substituting (3.9), (3.16) in (3.3), (3.4), we have for the external 
force system 


g \'? Zt-) 
ots | Ne. + | ( ) 


/2 


Zi- /2 
+7; (aS — (3.22) 


and for the external moment system 


1/2 
Me = f + ae 
h a 


ZY g /2 


| vo | (ay ° (3.23) 


It should be emphasized that the above expressions are exact, no approximation or 
assumption based on the thickness of the shell or plate being involved. 
Now we assume that E“ and X‘ can be expanded in power series in x°: 


«6m! 9 m! 


We also introduce the abbreviations 


= 3(4HaS + 58), C8 =Ka5+Hb8, DS = (3.25) 
d™ = — hy, ™ = hy, + d=d, =4#, (3.26) 
(-) 
Omi = ————- hts) ht), = 
| | we_,| (3.27) 
zi 
poi Tart ‘1+ tay Poi = pi, (3.28) 
| | | 


Here H and K are mean and total curvature, as in (3.13), and ¢ is the thickness of 
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the shell or plate, measured normally to So. We substitute (3.7), (3.24) into (3.18)- 


(3.23), and carry out the integrations. This leads to series in ¢, d, 


namely, 

= Etyt + (2HEto + Eft) + (2KEf + 4HEG + Ei) —- 

+ Ri (E*), 

T# = + + Efi) + (2C8 Ete + + Ei) 

+ Ri) (E*), 

4(3) 

Ler = + 2(BS Efe, + Eth) + RoE}, 


d 
= + (2HX{o + Xtu) [+ + 4H + Xie) 


+ + 2HQ‘? + + Ray(X®), 


(2) 


X fot + (BYXfo) + xt) > fo) + 


Fe 


+ Pe + + + + Ri) (X*), 


(2) 


Me 


+ + + C2Q@7 + 


(3) 


2(BEX + Xin) — Ris) (X*) 


In these series the remainders are as follows: 


t® 


Ra (E*) => 1 { KEim+2) 
m+5\(m+ 4)! (m+ 3)! (m+ 2)! 
(m+ 2)! (m+ 1)! 
Rw (E%) = rE By Etm+3) Efn+1) 
(m +3)! (m+ 2)! (m+ 1)! 
(m+ 3)! (m+ 2)! (m+ 1)! m! 
m+5 (m+ 2)! (m+ 1)! m! 
m+4\(m +2)! (m+ 1)! m! 
1 DS Eth+1) (m+6) 


+ 


(3. 


(3. 


(3. 


(3. 


(3. 


(3. 


(3. 


(3. 


(3. 


32) 


33) 


34) 


35) 


36) 


37) 
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M31) 
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These are accurate expressions, the summations being all for m=0, 2, 4,--+-. The 
remainders Ra)(X°), Ry (X*), Rys)(X*) are obtained simply by replacing Ej, with 
in Ray(E*°), with in Re (EZ**), Emm with X?, in respectively. 

In the ordinary case where the body force is the force of gravity, we can regard X‘ 
as a parallel vector field of constant magnitude. Then by (3.9), we have the following 
relations : 


Xfm) = 0 for m2 1, (3.38) 


Xf) = — Xfm) =O for m2 2. 


And consequently the tensor F‘, M¢ in (3.32)—(3.34) can be simplified to the (¢(, d“™) 
polynomial of few terms. The most important equations of this section, for future 
use, are (3.29)—(3.34) 

4. The tensors7°’, L“* in terms of the six quantities », 3, gas. We shall devote the 
present section to finding expressions for the tensors, (3.29)—(3.31), in terms of the 
six quantities pag and gas. Here pag represents the extension of the reference sur- 
face So and gag is closely connected to the change of curvature of the reference surface 
during the deformation; both were introduced in the paper [5] (p. 114, Eq. (44)), 
but will be formally defined below in (4.4). 

We shall now proceed to determine Ej?j, Eff, defined in (3.24) in terms of pag 
and gas. This is accomplished by means of (i) the equations of microscopic equilib- 
rium, (ii) the stress-strain relation, (iii) some geometrical relations, (iv) the conditions 
on the upper and lower boundary surfaces. The successive steps are as follows: 

(I) By means of (i), we express Ef) successively for m=1, 2,---+ in terms of 
Evo) and where n=0, 1, 2, - ++, (m—1). 

(II) With the aid of (ii), (iii) and the results of step (I), we determine Eg: Ett» 
Efi, EQ), Ej, successively in terms of Ef), Pas, Jas: 

(III) Then using (iv) and the results in step (II), we determine Ej, in terms 
of Pas, as. The surface force system (P', Q') is supposed to be given. Thus we have 
at once Ej§) and Eff for all m in terms of pas and gas. 

(IV) Substituting Eg, from step (III) into (3.30), (3.31), we obtain the required 
expressions for T** and L**. The expression for T*° can be found either by substi- 
tuting Eff) from step (III) into (3.29), or by using the equation of macroscopic 
equilibrium (2.5c). 

The geometry of strain and the definitions of pag and das. Let us introduce comoving 
coordinates [7, 8]. The same coordinates are attached to each particle during de- 
formation, and the coordinates form a normal system in the strained state. The funda- 
mental tensor in the strained state is g;; (satisfying (2.1), (2.2)), and in the unstrained 
state it is gj. Let Sj be the surface in the unstrained state which is carried over into 
the reference surface So in the strained state (after we reach section 6, we shall define 
Sj to be the middle surface in the unstrained state). The parametric lines of x° are 
not in general normal to the reference surface Sj in the unstrained state. 

The strain tensor e;; is defined as 


ei; = (gis — ii); (4.1) 
this is the definition usually adopted (cf. [7, 8, 9]). For small deformation, the prin- 


cipal part of the extension of an element in the direction A‘ is e=e;; A‘ 4. 
We shall throughout raise suffixes by means of g‘/; thus 


J q 
J 
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= gmignigt;, gm = giigim, (4. 2a) 
emn = gmignie; ;, = gle im. (4.2b) 
Now we assume that e;; is expansible in power series in x®, so that 
> — (4.3) 
m=0 
Let us define 

Pig = Clojiis Qij = Tig = (4.4) 


obviously these tensors are symmetric. The six components Pag, dag are the basic de- 


pendent variables of our theory. 

For small strain »;; must be small, yet qi;, ri;, etc. may be finite in a thin shell or 
plate. These quantities are not independent but are connected by certain geometrical 
relations. Since space is flat, g;; and gj, are not arbitrary functions of the coordinates. 
They must satisfy the equations 

= 0, (4. 5a) 
R: = 0. 
(4. 5b) 
Here is the curvature tensor for 


= 3( git, + — — 
[il, m],[jk, n], — [ik, m},[jl, (4.6) 


while Rj: is the corresponding curvature tensor for gj, 


= A(gin, ix + Sik, it 
+ [il, — [ik, m],-[jl, (4.7) 
It should be noted that R;jx: is the curvature tensor in space, while R,asy is the two 
dimensional curvature tensor of So (cf. (3.11)). In (4.7), g’"" by definition denotes the 
cofactor of g,,, in g’, divided by g’; namely, 
1 


— 
g = 


where e”*! is the usual permutation symbol, and g’ is the determinant of gj, 


g = (4.9) 
It follows that 


where 
= (4. 11) 


| 

4 

‘i 


1944] INTRINSIC THEORY OF SHELLS AND PLATES 307 
this is a contravariant tensor and satisfies the relation 
= 5154 — 64,5}. (4.12) 


The equations (4.5a, b) form a set of twelve independent geometrical conditions 
for gi; and g;,. We shall now regard g;; as given, and develop (4.5b) explicitly in terms 
of e;; by using (4.5a), (4.1). The resulting expressions represent the six independent 
geometrical conditions for e;;: 


+ 2e""{ [il, m],[jk, n], — [ik, m],[jl, n],} 
— 2gmn{[il, m],[jk, mle + [il, m].[7k, — m].[jl, n], — m],[jl, n].} 
+ + + = 0. (4.13) 
This is a polynomial in e;;, with linear terms explicitly stated. The other terms pro- 
ceed with increasing degree in e;;, and have the following exact expressions: 
+ [il, m].[jk, n]. — [ik, m].[jl, n].} 
+ A(emer — eet) { [il, m],[jk, n], — [ik, m],[jl, n]o} 
+ — em) { [il, m].[jk, nm], + [il, m] n]. 
— [ik, m].[jl, n], — [ik, m],[jl, n].}, (4. 13a) 
— pw( Git, ik,it — Bik,it — Bits) 
— — { [il, m].[jk, n]. — [ik, m].[jl, n].} 
— [il, m] .[jk, — [ik, m].[jl, n], 


+ [il, m],[jk, n]. — [ik, m],[jl, n].}, (4. 13b) 

= pw( Cit, sk + — jt — (4.13c) 
[ij, = + — Bein) (4. 13d) 
[ij, Rle = + — (4. 13e) 


Three of (4.13) involve only e;;, ei; and their derivatives with respect to x*. 
When x°=0, these three conditions become the equations of compatibility of the 
twelve quantities ~;; and q;;; a detailed formulation of these equations will be given 
in section 5. 

The other three of (4.13) involve not only e;;, €i;, and their derivatives with re- 
spect to x*, but also éas,00. In fact, by these equations, we express €4,00 in terms of e43, 
é:j,0 and their derivatives with respect to x*. When x°=0, these equations give us the 
expression of rag in terms of p;;, qi; and their derivatives with respect to x*. Putting 
i=a, j7=0, 1=0 and x°=0 in (4.13), and solving for we obtain 
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Yas = = + — Poojas + + baz) graze + 
— + + — bhaz — prow 
— Jo0Jas — + PG), (4.14) 


where O,)as is a residual term, to be explained below. 

To find e;3}2s, we have to differentiate @ag,o0 with respect to x°, and put x°=0 in 
the resulting equations. By (4.14), this gives e;3)2g in terms of jj, gi;, rio and their 
derivatives with respect to x*. To find the other eéjmjas, we merely repeat the process 
over and over again. Thus we can express @jmjag in terms of pij, Qij, €{njco and their de- 
rivatives with respect to x*, where n <m. 

Actually, to obtain the principal parts of the final equations, we only require 
explicit calculation for m=2, as in (4.14). 

All the above results are of a purely geometrical nature. 

The symbol O.2)(p?, #@) in (4.14) represents the terms which are not explicitly 
calculated. The quantities in parentheses indicate order of magnitude of these terms 
for small # and g, which denote the magnitudes of the tensors pi;, gi; respectively. 
If #, g approach zero simultaneously and independently, these terms converge to zero 
at least as fast as $g or f*. Symbolically, we may write 


£G) = Oas(p?) + (4. 14a) 


The label “(2)” is to distinguish this O-symbol from later symbols of the same type. 
The indices attached to Oy) are the tensorial indices of every terms involved. We 
assume throughout that differentiation with respect to the coordinates x* does not 
change the order of magnitude; i.e., pas, Pas,y are of the same order. On the other 
hand, we never make any assumption regarding the effect of differentiation with re- 
spect to x°. 

This O-symbol notation will be used extensively throughout the paper. The nota- 
tions used inside the parentheses of O-symbols are given in the following table: 


Tensors | bag | bag | Cop | | | ES | Pii Pap | Pap | | 
Symbols for magnitudes b | b e | é| E | EB >|? | Pp | P| Pp | 

Tensors ap | Qi | Xie | 
Symbols for magnitudes q | q 0 | Q | _ | X 


The expressions of Ejp, in terms of Ejj, and Epi (n=0, 1,-+-, (m—1)). 
We start the process outlined in the beginning of this section by writing down the 
microscopic equations of equilibrium under the body force X?: 

E*iy, + Xi = 0, (4.15) 
the double stroke indicating covariant differentiation using g;,. Putting in turn 7=0 
and 1=a, we get three equations which may be written as 

= — 3g — — X°, (4. 16) 
g 


3 
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where g under the stroke indicates covariant differentiation using gas, 


= + BEY, 


(4.18) 
a 
Et¢,, = E™*, + { + *, (4.19) 
9 Aw) @ Ar) 


This operation should be clearly distinguished from that indicated in (2.6). Putting 
x°=0 in (4.16), (4.17), we get 


Ett; = — — + ER, — (4. 19a) 
Ef}, = — — — — (4.19b) 


These are the expressions of Efi, in terms of Ej) and Ej, the body forces being sup- 
posed given. To find Ej), we have to differentiate (4.16), (4.17) with respect to x°® 
over and over again, and put x°=0 in each of the resulting equations. We shall only 
write the expressions of Ej} explicitly: 


Et) = — Xt + 2HX to + + + — anKEM 

+ (8H? — + bX + 4(H p, (4. 20a) 
Ets, = — Xfi + (0% + fy — + (OF + Ein 

— 3(4Hat + ER + (10Hb¢ + — 8Ka%) (4.20b) 


This completes Step (I). 
Stress-strain relations and expressions for Ejm, in terms of Pas, Qas, 
We shall accept, as the stress-strain relation for an isotropic body, the equation 


1 
(1+ — 20 


ij 


5 { + (1 — + O'%,(@), (4.21) 


where @ is Poisson’s ratio, and E* is of course the reduced stress (see section 2) ; there- 
fore Young’s modulus does not appear. For small strain problems, @ is small, and 
the terms in O*/(é?) are negligible in comparison with the terms linear in e* in (4.21); 
in that case, (4.21) becomes the usual linear stress-strain relation for an isotropic 
body (in rectangular Cartesians, see [1], p. 102, Eq. (18)). Any modification, such as 
replacement of g‘’ by g’*/, leads to no real change, because the difference is taken 
care of by the O-symbol. 

From (4:21), we have 


1 
(1 + ¢)(1 — 20) {og + (1 — c)eoo} + (4.22) 
1 
(1 + ¢)(1 — 2c) tog®eoo + + (1 — 20) g**g Jee} + (4.23) 
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It is evident from (4.21) that all the stress components are small, of at least as high 
an order as the largest strain components. The elimination of eo; from the last three 
equations gives 


1 
= — + (1 — Jen + Oe, Be). (4.25) 
= 


When x°=0, the equation (4.25) becomes 


o 


[0] At) Par + (p*, Eop), 


E%, = 
(4.26) 


where the abbreviation Ag means 


1 
Att = feartar + (1 — (4.27 


Equation (4.26) is the required expression for Eva: 
The substitution of Ea from (4.26) into (4.19a, b) gives 


2(1 — 20) 
(1 — o) 
0%, (bp?, bEop, bF3), (4.28) 
Eft} = — Xfo — — — a** — 
+0%8,(p?, pEo, E3). (4.29) 


Here Xj) are supposed to be given. These are the required expressions of Ej). 
Now, differentiating (4.25) with respect to x® and putting x°=0 in the resulting 
equations, we obtain in consequence of (4.26), (4.28) 


8 ‘ 2(1 — 20) 
= — a®X%) — 4 
+ a8, = ary __ a™®A 
2 — 0) (1) (1) wy 
+ AGPgn + (4.30) 
where 
= (bp, bE?, Pq, pX, XE,). (4.31) 


Equation (4.30) is the required expression of Eff. 
The substitution of EG: Ei from (4.26), (4.30) into (4.20a, b) gives the required 


expressions for E73}: 
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2(1 — 20) a 
Eth) = ———— — Xtu + 
(1 — 
2(1 — 
+ {2(2 — 3c)H? — (1 — o) K} 4 
2(2 — 3e) 
— 30) 
(4 — 30)Hb™ + Ha" — (3 — «)Ka™ $ pn, 
1 — o? 
1 
+ ———— [{4oHa™ + (1 — o)b™} + 0%). (4.32) 


2(1 — oa?) 


20 
a** (HE +- (10Hb¢ 8 Kaz) E fo) Aa)" Qed 8 
—o 
2c? 
— 2 
4+ + 4 Ha®a™ 
(1 — o) 
+ (1 — + + pei + (4.33) 
Here 
8) = pEo, Ei, bpg, bgEo, bpX, LXE»), (4.34) 
&) = bpEo, bE, pq, qEo, pX, XE»). (4.35) 


To find Ee: we have to determine pio, gio and then rag in terms of pag, as, Ey 
Eliminating éo9 from (4.22), (4.24) and solving the resulting expression for é@ao, we 


obtain a 
Cad = (1 + 7) + O10)a(e?, eE, E?). (4.36) 


when x°=0, this becomes 
pao = (1 + + pEo, E2). (4.37) 


We now differentiate (4.36) with respect to x° and put x°=0. In consequence of (4.26), 
(4.28), (4.29), this gives 


| 
4 
q 
| 
ia 
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i+o 
{(1 — o)Xpoja + CEM + 21 — 0) — 


(4.38) 


— 


@* + 0 11) Fi, Pq, pX, XE,). 
Similarly, from (4.22), the values of éo9 and éo0,9 on x°=0 are, by (4.28), (4.29) 
(1 + — 2¢) 


Poo = Et) + pEo, 3), (4.39) 
(i= 4) 1— 
(1 + «)(1 — 2c) — 29 
(1 ) l-—o i- 
1 
{40(1 — 20)Ha™ + (1 — o)b™} per + 
where 
O.1s)00 = bpEo, bE’, pq, qEo, pX, XE1). (4.41) 


Substituting pio, gio from (4.37)—(4.40) into (4.41), we obtain 


rap = 1B (q) + Bp) + Q(X) + 1B (Eo) + 19g) + Oanas, (4.42) 


where the abbreviations represent 


1 
8 (9) 8 + bga ) (4.428) 
(1 + 
ros (X) = — {2(1 — + + (1 — 20)bapX (4.42) 
l+o (1 — 20)? 
(1 — a) 


— 2(1 + + (HEF) 10 + bias, 


(1 — 2c) 
— — + (4.42d) 


r3 (¢) = ; — que (4.42e) 
= pEo, Ei, Pq; qEo, PX, (4.42f) 


It will be noted that r{4(q) is linear in the q’s, r2}(p) is linear in the p’s, and so on. 
We now differentiate (4.25) with respect to x° twice, and put x°=0 in the resulting 
equation. In consequence of (4.28), (4.32), (4.42), this gives 


| 
1 
| 
| 

| 

‘ 
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Es, = + (p) + + + + Os, (4.43) 


Otis ork, (p’, pEo, Ei, Pq gEo, px, XE,, qX); (4.43a) 


the other abbreviations are 


1 
E%},;* (q) — { + Saba + 3(1 — o)(b**a + Sans, (4.43b) 
(p) = — + 
+ (1 — + a%ahr)a™} 
40? o(7 — 9a) 
1 — a? (1 — o) 
o(7 — 120) (1 — 
 Hat1h8 + ————_ — 3(1 — o) 
(1 4 
+ 6(1 — 0) + — (1 — (4.43) 
1 
— at + (4. 43d) 
1 
= — { (180? — 140 + 1)Hb* 
— — + 20(4 — 50)Ka%®} E%, 
1 
1 
+ + ab, — — (4.43e) 


Equation (4.43) is the required expression of Ef: 

By similar steps, we can express all the other Ej and Ee in terms of pag, Gas 
and Ejs. However, for the purpose of the first approximation in the equations of 
equilibrium, the knowledge of the required expressions of Ej), Eji;, Ej is sufficient. 
This completes Step (II). 

Conditions of the boundary surfaces, and the expressions of Ejo, in terms of pap 
and qas. We have now succeeded in expressing Ef, in terms of pap, Jas, Ej. Our next 
task is to express Ej in terms of pas, gas by applying the following six surface con- 
ditions: 


where | 
4 
4 
; 
| 
| 
r 
i 
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for = + ha, (4.44a) 
= Zi, for = — hy». (4.44b) 


The unit normal vectors N;4):, Nii, drawn out from the shell or plate, are deter- 
mined as functions of x* by the equations 


N N => N a (4.45a) 
Neyo = (1 — » = (1 , (4.45b) 
the positive roots being understood. Z,4):, Z;-); are tensor components of the given 


loads per unit area applied to the upper and lower surfaces. 
Making use of (4.45a, b), we find that (4.44a, b) take the form 


| | 


Substituting (3.24) into (4.46), and adding the resulting expressions, we obtain in 
consequence of (3.26), (3.28) 


EQ) — = + = (4.46) 


+ Eid + — Bt 4... 
1 


1 
~ Ezz + + — | (4.47b) 


Substituting Ej; from the results of Step (II) into (4.47a, b), we obtain a set of three 
equations in terms of nine quantities pas, gas, Ejq. Solving these equations for Ej, 
we have 


1 — 20 
an {(2 — 3c)2H? — (1 — o)K}Q%® — |, 


2(1 — 20) 
1 
20 — 0) { 20Ha™ +(i- o)b*} (pay 


— 3c) 


— (3 o)Ka™ pat 
1 — a) 


1 
4(1 — o?) 


8(1 {4Hoa™ + (1 = o)b™ } + Ofte)» (4.48a) 


KS 
4 
4 
( 


i944] INTRINSIC THEORY OF SHELLS AND PLATES 315 


= + + 3(6% + 


4(1 — o) 
+ 3X fod — Toit) 
{ (4Ha® + — + 2Ha®)Q?, } 


+ 3(SHbg + 4H%af — + 3{ — + 


1 {| (b%#8 + } hh) 20? H aB 


(1 — 
(b=* oi) | pase 


— + + (4.48b) 


Here Offs), Off, stand as usual for terms not explicitly calculated. It is possible to 
exhibit their orders of magnitude as in (4.43a), for example, but the expressions are 
very long and will therefore be omitted here; the full expressions may be found in the 
author's Ph.D. Thesis, “The intrinsic theory of elastic shells and plates” (University 
of Toronto Library). The important fact about these residual terms is that, in the 
case of small strain, they are small compared with the terms shown explicitly. 

The substitution of Ej, from (4.48a, b) into (4.26), (4.28), (4.29), (4.30), (4.32), 


(4.33), (4.43) etc., gives the expressions for Ej) in terms of pag and gas. This com- 


pletes Step (III). 

of T*°, in terms of pas, Jas. If we substitute the expressions 
of Eff, from the veal of Step (III) into (3.30), (3. 31), we immediately obtain the 
set of eight quantities 7** and L** in terms of the six quantities pg and gag; the quan- 
tities Xj), Q' are supposed to be given. Therefore for the membrane stress tensor 
T** and the bending moment tensor L**, we have 


T# = Athy part BOS + + 24 


+ CO" Patt — (aA + aA) + 


o 
2(1 — o) + TE TH {X} TS {Q} + Offs), (4.49) 
L® = prsd® + BE prst® + 44 grat } 


+ a5(30°d — 4X jot — 


+ Oth), (4.50) 


# 
q 


316 WEI-ZANG CHIEN [Vol. I, No. 4 


where Offs), Offs) stand for the residual terms. The other abbreviations are 


1 
= + (1 — 
(1) ( (4. 51a) 
AQ) = ca + (1 — — 
6(1 — o?) 1 — o? 
Ba; = {4o0Ha*™ + (1 — o)b™} (4.51c) 
8(1 — o)(1 — o?) 
BG) = ———— {o(8Ha*a" — — 
12(1 — o?) 
+ (1 — — — (4.51d) 
Bay = ———— — — 
4(1 — o?) 
+ (1 — — — (4. 51e) | 
Bo = Hat*a® + (1 — 0)4Ha®a™ — + 
6(1 — o?) 
— 2(1 — o)(b®a"™ + (4.51f) 
Cay. = —— 4 (4 — 30) 4+ ————— H?a™ 
4(1 — o?) (1 | 
@ cas 1 { 21 + argh) 1— sen 
6(1 — o?) 
4o? a(1 — 4c) 
— $(1 — o) + H?a®%a” — + ————— Hb**a™ 
a(1 — 3c) 
— 8(1 — — Kastan} (4. 51h) 
2(1 — 2c) 1 
Tt {X°} = ——— ont — —— Axi a® X 


1 
61 {4(1 40) fo) + 4o0Ha® ote (4.51i) 


4(1 — a) 6(1—a 
— P(X + X (4.51}) 


> 

| 
A 
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8(1 4(1 — a) 


(i 


1 


40 — 1 12¢ 
p28 
l-—o 


+ (4. 51k) 


Furthermore, by solving (2.5c), we have 


= + (4.52) 


in which L* is given by (4.50) and M® by (3.34). 

Equations (4.49), (4.50) express T** and L*8 in terms of pag and gas. When L** 
is known, T°° is calculated from (4.52). This completes the last step of the procedure 
outlined at the beginning of this section. 

It should be noted that pag and gag correspond respectively to the extension and 
change of curvature of the reference surface So; Xfm is the normal derivatives of the 
mth order on Sp of the body force, supposed to be given. If the form of the reference 
surface So is given (in the strained state), the following quantities are known: @ag, 
(1/2)bas, (1/2)¢ag are the first, second and third fundamental tensors, H the mean 
curvature and K the total curvature. Furthermore, if we know the positions of the 
boundary surfaces of the shell in the strained state and also the surface loads on the 
boundary surfaces, then the quantities ¢™, d™, Q™‘, P(/™# can be determined by 
using (3.26)—(3.28). 

5. Equations of equilibrium and compatibility in terms of the six unknowns pas 
and gas. Having now expressed the macroscopic stress tensors in terms of pas and Gas 
we shall substitute these expressions into the macroscopic equations of equilibrium, 
(2.8a, b). In consequence of (3.32)—(3.34), this gives the following three equations 
in terms of the six unknowns pag and dag: 


— part + AAT (Geral) + 
— Bs get + — bp 
— 3b,,Buy prd® + | xt 
pat + FAD (pard™) + (BR + 


HOt + + KP® +4 7%, {Or} 4+ 1%) {Pr} + 1%s){ 


+ {X*} = (5.1) 


a 
it 
o 
20 
is 

| 

1 
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a ra Ta 
+ (grat) + (Boy pryd™) 1, + Prt) |, 


2 a Ar ad 


+ 1p — 1p + 


4+ — (00%) + + — + 1%) {QF} 


2(1 — o) 
+ I~2){X°} + = Ob. (5.2) 
The symbols Ag”, Ag™, given by (4.51a, b), are functions of das; ne. ne, 


Bee, B®” given by (4.51c, d, e, f), are linear functions of bag; and Cah", Ci, 
(4) (5) » 8 y B (1) (2) 


given by (4.51g, h), are quadratic functions of bas. The other abbreviations in (5.1), 
(5.2) are listed as follows: 


o 
1°, {Q7} = ———— A { (BF + 2H } —- ———— HQ" det 
{ (oF + | 20 — Ord, 
2(1 — o) 6(1 — 
o 40 — 1 
6 1-—o 
+ (HD) + {QM + (Kak + (5. 3a) 
Pt} = {3(4Hax + bX) P + KP} (5.3b) 
1%) {X°} = Xfoyt + (2HXto) + — — rz 
2 3(1 — o) 
1 
{ (4H? + 60K) + 2(1 — o)HXty + fe 
q 
2(1 — 2c) 
2(1 — o) l-—o 


{ Xt} = + + + Je 


“di 


Ti {Q"} = + 2Hax + 207d) xt 


o 
l—o 


: 

q ‘ 
a 
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20 


1 


40 — 1 120 


| 
18 


+ (2Haz + + (Kb; + 2H; — 2KHaz)Q*, (5. 4a) 


1%) = + 4Xfyd™ — + tod + 


+ — { (Xpoyen + } ip 
+ 3(2Hat + (5.4c) 


Equations (5.1), (5.2) are the three equations of equilibrium in terms of the six un- 
knowns Pag and gag. All the other quantities are supposed to be given. 

The basic unknowns are the six quantities pag and gag. To find them, we have to 
solve the three equations of equilibrium (5.1), (5.2), together with the three equations 
of compatibility, funished by the geometrical conditions. 

We shall now convert (4.13) into a form in which pag and gag are the only un- 
knowns. Let us put in turn i=0, j=a, k=8,1=y, andi=p, j=a, l=vy in (4.13), 
and then put x°=0; we thus obtain three equations in p;; and g;;, We now substitute 
(4.37)—(4.40) for poo, Pas Goo, Yao, and (4.48a, b) for Eq; this gives three equations in 
Pag and gas. Two of these, after being multiplied by nj and simplified by using (3.13), 
(3.14), can be written as follows: 


+ (Pariy + Prria — 


|} = (5.5a) 


nio) {24000 


+ (1+ 0) [Ona — GapKQy + 


The third equation, after being multiplied by nfjnf}, becomes 
2(1 — 


Posiay — 110) 210) Jay + 


2 — 
(1 + «)(1 — 2c) 
1—o 


OK — (1 + o)(b** = 23). (5.5b) 


1 
} 

4 

§ 
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Equations (5.5a, b) contain the three equations of compatibility in pag and gag; these 
are to be associated with the three equations of equilibrium (5.1), (5.2) for the solu- 
tion of plate and shell problems. 

6. The equations of equilibrium and compatibility referred to the middle surface 
in the unstrained state. The quantities aus, bag, which occur in (5.1), (5.2), (5.5a, b), 
refer to the strained state. But it is usual in elasticity to regard the unstrained state 
as given, rather than the strained state, and from this point of view @as, bag are un- 
known. We should use, instead of them, the fundamental tensors of that surface Sj 
in the unstrained state, which passes over into the reference surface So in the strained 
state. So far So has been quite general. Now we shall follow the usual method by 
choosing So so that Sj is the middle surface of the shell or plate in the unstrained 
state. This means that Sp is not accurately the middle surface in the strained state. 

For the metric in the unstrained state, we have by (4.1) 


= Bij — 20:33 (6.1) 
and so, if we define 
we have 
Gas = das — 2pap, = — 2pao, = 1 — 2poo, 
(6.3) 
Das bas = 24 a8; bao = — 2¢a0; boo = — 2q00, 


where /;; and qi; are defined as in (4.4). aig is the fundamental tensor of Sg ; we re- 
gard it as given. We may substitute in the proceding theory 
ap = ap + 2pas. (6.4) 


The tensor bi, does not represent the curvature of the middle surface Sy in the natu- 
ral state; since in general g’,~0 on Sj, the parametric lines of x® cuts Sg obliquely. 
Let us introduce the normal coordinates x‘ based on Sj , choosing x*=x* on Sg, and 
x° normal to Sj. The metric g;; corresponding to the coordinates x’ satisfies 


&o0 = Bao = 0. (6.5) 
Let us put 


98 « 
Aap = (Bap) bas = ( (6.6) 


Now 4:; and gi, are metrics corresponding to the coordinate systems x‘ and x‘ re- 
spectively, for the description of the geometry of the unstrained state. Hence the 
tensor g;; determines g;,, and vice versa. 

It should be noted that the quantities (1/2)b.ag are the coefficients of the second 
fundamental form of Sj ; they vanish if Sg is a plane. The radius of curvature R in 
the direction of a unit vector wf (R counted positive when Sj is convex in the sense 
of x° increasing) is given by 


bas (6.7) 


4 

4 

: 

4 

R 2 
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We proceed to find in particular the quantities (6.6). In the first place, since 
x*=x% on Sj, we have 


aap = a8; Gap = Aap + 2pas (6.8) 


We now follow a straight line normal to Sj, starting at the point x*. Since it is a 
geodesic, we have (cf. 6, p. 301) 


1 dxi dx* 
+ { 0, (6.9) 
(dx®)? dx® dx® 


and so we can develop x‘ as power series in x®°. For x°=0, we have 


dxi a 
= x*, = — 


dx®  (g/%) 1/2 (6.10) 


where 2’! is the conjugate of gj. The last follows from the fact that the line is nor- 
mal to Sg, and so 
dx! dx’ dx‘ 


bai = 0, gig — = 1. (6.11) 


On carrying out the development in power series, and using the transformation 


435 = fa. (6.12) 
we obtain after a little calculation 


bas — 2gas = = + atop + + (6.13) 


where @’‘i=Z'* for x°=0. The Christoffel symbol is calculated for aig. By (4.10) and 
(6.3), we have 


a’ = 1 + 2poo + = 2p + (6.14) 
Thus (6.13) becomes 


bas bis bag baspoo 2Poals 2Posia + O (26) a8( (6.15) 


where a under the stroke indicates the covariant differentiation with respect to x* 
and @ag. 
Let us define pag, as SO that 


2q = bap — bas, 2Pa8 = 2pas = dap — Aas; (6.16) 


then the extension and change of curvature of the middle surface Sj along the direc- 
tion of a unit vector pio) are given by 


Pas (6.17) 
From (6.15), (6.16), we have in consequence of (4.37), (4.39), (4.48a, b) and (6.8) 


q 
4 
i 
a a 
as 
j 7 
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(1 + — 2c) 
= 2qas — + bagQ® + (1 + o)Qais 
1-—o 2(1 — o) a 
+ (1 + + O(27) 08, (6.18) 
bas = 2qas + bas, Gap = 2Pap+ aap, pap = Pas. (6.19) 


The symbol O72 represents the residual terms with the order of magnitude shown 
symbolically by 

O (27) a8 O (27)as(p”, Op, pd, bO%®, bQd, Qt, qQd, 

bX°d, X%, Xd, Xt, qXd). (6.20) 


Let us now denote the thickness of the shell or plate in the natural state by 2h. 
Then by (6.10), (4.10), (6.1), we have 


00 0 2 0 
h= Genin -f [1—e +O s)(é@)]dx , (6. 21a) 
dx® 00 00 0 
h -{ -{ 1— O 6.21b 
By (4.3), these become 
00 
h = heyy — hey +O 9) herr, (6.22a) 
00 
h=hey— p + GPA). (6.22b) 


Substituting (4.37)—(4.40) for pio, gio, (4.48a, b) for Ejg, and (6.18), (6.19) for gas, Pas, 
Gas, bag into (6.22a, b), we have two equations for the determination of h,,) and h,_). 
We now solve these equations for h,,) and A). In deciding what terms to retain ex- 
plicitly, we note that, for a thin or plate undergoing small strain, the quantities h, 
Pas, Q', are small. We obtain 


(1 + — 20) 


h =h- + 6.23 
(+) + 2(1 — Oh + ( a) 
(1 + — 2c) 
hy) = h = a™pwh Oh 0:32). (6.23b) 

2(1 — a) 


The O-symbols represent the residual terms with orders of magnitude shown symboli- 
cally by 


p*h, O°h, pOh, Xh?, gh’. (6.24) 

Hence from (6.23a, b), we have immediately, in the notation of (3.26), 

no + — 20 
im = 2(1 + ( ) + 0433), (6. 25a) 
1—o 2(1 — o) 
d™ = Or, (6.25b) 
where 

= Oss)(p*h", pOh", Oh", (6. 26a) 
= Ow (Xh™, gh", ph**). (6. 26b) 


J 
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Similarly, substituting (6.23a, b) into (3.27), (3.28), we get 


= + (6.27a) 
= Oh" + (6.27b) 
where 
Oiss) = Oiss)(Pph", POR", Pgh"), (6. 28a) 
Oiss) = Ofs)Oph", Ph", OXh", Ogh"*), (6. 28b) 


where P‘ and Q‘ defined by (3.28), represent the sum and difference of the components 
of the surface loads with respect to the surface So in the strained state. 

With (6.18), (6.19), (6.25), (6.27) established, the expression for T**, L«8, T°° in 
(4.49), (4.50), (4.52) will now be reduced to forms involving pas, Gas instead of pag, 
Gag. The results are as follows: 

(i) The membrane stress tensor T°, 


(ii) The bending moment tensor L*°, 


40 
+ fa HQ — — + 0%). (6.30) 
6(1 — o) 1-—o 


(iii) The shearing stress tensor T°, 
= 2{ AUS)” + 1 + 
a 


+ 3(4HP* + + + 


40 


— — + 4{(4Hat + + Xf 


+ + a% qn) + Of. (6.31) 
The residual terms in (6.29)—(6.31) are 
= O%s)(bp*h’, pOh*, XOh*, gph’, Xph’, bph*, gh’, 
Oh’, Xh'), (6. 32a) 
= 0%.,(p*h, Oth, Oph, ph’, bOh*, qQh', baph’, q*ph’, gh’), (6.32b) 
= O%%(Oph, OOh, OXh*, Ogh', bPph?, qPph’, bPQh*, pgh, 
bPgh!, Pqth®, bPh', bph’, gh’, Xh'). (6.32c) 


j 
4 
1 
om 
4, 
= 
aa 
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The abbreviations are as follows: 


aBrr 1 
Ajay = 1 {ca%®a™ + (1 — c)arta®}, (6.33a) 
— 1) 
AG’ = ——— + (6.33b) 
3(1 — a) 
Ag = — ath, — — (6.330) 


amrdpy 
A (5) 


o 
— 2a — + (6.33d) 


nfo) = a=det. (acs), =? =0, (6. 33e) 
= — 1, 
H = ha™bn, (6.33f) 


where H is the mean curvature of the middle surface Sg in the unstrained state. All 
these quantities are determined by the geometry of the middle surface (in the un- 
strained state) of the shell or plate. 

With (6.18), (6.19), (6.25), (6.27) established, (5.1), (5.2) can be reduced by direct 
substitution to the form involving pag and qag instead of pag and gas. Thus we have 
the following three equations of equilibrium in terms of pag and as: 


a 


— + AC + AG 


P 2(1 — 2c) 
+ + + (OPA) + 
a 
1 — 2¢ 
+ 1 = Oo), (6.34) 


(per) 1p + AQ |p + is 
A a a 


— (qryqnsh®) |p + (Qrsh®) + a%®(h) 
a a a 


+ P* + 2X + (2Hat + 
+ (a®qnaz + = Of), (6.35) 
where 
0%) = O's0)(Oph, 00h, OXh?, Ogh?, bPh?, qPh?, q*?P*h’, qbP*h?, 
bph*, b0°h', Oh’, Xh', gph, gh’), (6. 36a) 
0%1) = 0%1)(p*h, Oph, b*Ph?, bqPh?, q*Ph?, ph’, Xh', qphi, gh®). (6.36b) 


: 
— 
es 
; 
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The abbreviations Agr’, are given by (6.33a, b, c), while AQy™ is 
given by 


a’ta™a™ — (2 — 


1 — 4)o 
A’e) 
6(1 — o?) 


+ 9(1 — + — (6.37) 


All these tensors are determined by the geometry of the unstrained state. 

If we were to substitute (6.18), (6.19), (6.25), (6.27) into (5.5a, b), we would im- 
mediately obtain three equations of compatibility for pas, das, With certain terms not 
explicitly calculated. However it is wiser to adopt an entirely different method, for 
the required equations of compatibility can be obtained in an exact form by a purely 
geometrical method. _ 

We first introduce the equations of Codazzi and Gauss for the reference surface 
So in the strained state, 


bapiy — = 9, (6. 38a) 
= — beybpa), (6.38b) 
and the corresponding equations for the middle surface Sj in the unstrained state, 
basiy — Baris = 9, (6. 39a) 
a 
Ryasy = — byybpa). (6.39b) 
Here we recall that 
bapiy = bap,y — a*?[ary, plabse — pladar, (6.40a) 
Dasiy = — a*[ary, — (6.40b) 
+ af [oy, — [08, Ja}, (6.40c) 
= + Aas.py — — 
+ a™{ [oy, r]a[aB, — [08, relay, dJa}. (6. 40d) 
Furthermore by definition, we have 
at = = ay (6.41a) 
a™ = a= (6.41b) 
ay = by = fee" ay = = (6.41c) 


where 6; is the Kronecker delta. Substitution of a,s from (6.19) into (6.41a, c) gives, 
with (6.41b) 


a 

av’ = (2nfointapas + a), (6.42a) 
a 

1 + + 2pna™, (6.42b) 
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a 
= — [85 + + + 
where nj) is given as in (6.33f). 
We now multiply (6.38a), (6.38b) respectively by (a/a)nf}, (a/a)nfiinfj, and sub- 
stitute das, bas from (6.19) into the resulting equations. This gives, in consequence 
of (6.42a, b, c) and (6.39a, b), the following three equations in pag and qas: 


a 


nid) + bsx)(a®™ + nfs) (Pauly + Pywla Pay|w) = 0, (6.43) 


(1 + prspey + 2a pry) { nfo) 
a 
+ it} dow + pasK — — b*)qas} 
+ (a™ + Pas) + — Proit)(Pawis + Paola — Papi) = 0. (6.44) 
a a a a a a 


Equations (6.43), (6.44) are the three equations of compatibility for pag and qas; these 
are to be associated with the three equations of equilibrium, (6.34), (6.35) for the solu- 
tions of plate or shell problems. We note that K is the total curvature of the middle 
surface Sj in the unstrained state, and satisfies 


K= by; — b™b,). (6.45) 


Conclusion. Equations (6.34), (6.35), (6.43), (6.44) are the final forms of six 
differential equations in the six unknowns pag and qas. Here pag and as, as indicated 
in (6.17), represent the extension and change of curvature of Sj, the middle surface 
in the unstrained state. With pag and qag known, the macroscopic tensors T**, T°, 
L** can be calculated from (6.29)—(6.31). 

We recall that the tensors aag and (1/2)bag are respectively the first and second 
fundamental tensors of the middle surface Sf in the unstrained state; H and K are 
the mean and total curvature of Sj as in (6.33e), (6.45); @ is Poisson’s ratio; 2h is 
the thickness of the shell or plate in the unstrained state; P‘ and Q‘ represent the sum 
and difference of the surface forces on the upper and lower surfaces of the shell or 
plate in the strained state as in (5.28); Xj, are the normal derivatives of the body 
force on the reference surface Sy in the strained state as in (3.24). All these quantities 
may be regard as given. The covariant differentiations are calculated for azg and x*. 

We also note that the six equations (6.34), (6.35), (6.43), (6.44) are exact, in the 
sense that no terms have been omitted, but of course the residual terms, represented 
by O-symbols, have not been calculated explicitly. However, it will be shown in Parts 
II and III that in all cases of small thickness and small strain, the residual terms are 
small compared with those shown explicitly, and it is legitimate to neglect them in a 
first approximation. 


(To be continued) 
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ELECTROMAGNETIC WAVES IN A BENT PIPE OF 
RECTANGULAR CROSS SECTION* 


BY 


KARLEM RIESS 
Tulane University 


The analysis of electromagnetic wave propagation in a bent pipe of rectangular 
cross section, (x =R, x=R+a, y=0, y=d), is based on the Maxwell field equations, 
expressed in cylindrical coordinates (r, 6, y) (Fig. 1). As in the case of the straight 

- pipe,' the time variation is given by the exponential e’*‘, where w is the angular fre- 
quency. The angular variation is given by e~’, where = is the propagation constant 
for the bent portion. The equations may be written 


— DE, — + rjwpH, = 0 
0E,/dy — dE,/dr + jopHs = 0 
+ Ey + ZTE, + rjwpH, = 0, 
0 


— 2H, — rdHs/dy — rjweE, = 0, (1) 

dH,/dy — dH,/dr — jweEs = 0, 

rdH»/dr + He + 2H, — rjweE, = 0, 


r0H,/dr + H, — 2H, + rdH,/dy = 0, 
rdE,/dr + E, — TE, + rdE,/dy = 0. 


In (1) He, H,, Hy, Es, E, and Ey are the components of magnetic and electric field, 
€ is the electric inductive capacity, and yw the magnetic inductive capacity. The elec- 
trical conductivity, o, and charge density, p, are assumed to be zero. 

The field components H,, H,, E, and E, may be expressed in terms of Hy and Es by 
various combinations of the equations (1). These give 


H,(Gr?) = — 2rdH»/dr — THs — r?jwedEs/dy, (2a) 
H,(Gr?) = — + rjweEs + 1? jwed Es /dr, (2b) 
E,(Gr?) = — ZrdE,/dr — + (2c) 
E, (Gr?) = — 2rdEs/dy — rjwpHs — (2d) 
where 
(Gr?) = + 


Using the last two of Eqs. (1) and Egs. (2), H,, H,, E, and E, may be eliminated 
and equations in H, and £, readily obtained. 


=| 0(rHe) rjwe OE, 
OrLGr or Oy? 


= 0, (3a) 


* Received June 29, 1943. The paper constitutes part of a thesis submitted in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy at Brown University, May, 1943. The author 
wishes to thank Professor L. Brillouin of Brown University for his many helpful suggestions. 

ai 1 This case has been discussed by Lord Rayleigh, Phil. Mag. 43, 125-132 (1897); Brillouin, Rev. Gén. 
’ de I’Elec. 40, 227-239 (1936); Schelkunoff, Proc. Inst. Radio Eng. 25, 1457-1492 (1937); Chu and Barrow, 
pe, q Proc. Inst. Radio Eng. 26, 1520-1555 (1938); Slater, Microwave transmission, McGraw-Hill, 1942, 
pp. 124-150. 
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odO(rE rj 0H, dG 
OrLGr or dy? = oy oF 
The boundary conditions for this case are 
y=0, y=b. Ee=0, £,=0. (4) 
r=R+ a4. E,=0, E,= 0. 
A more useful form may be obtained from (2d): 
y = 0, E, = 0, E, = 0. 
(5) 


7 = r=R+a. Ey = 0, O(rH¢)/dr = 0. 

By considering Ey and Hp as functions of y and r, these conditions establish the 
dependence of E, on sin kyy, and Hy on cos kyy, where ky =n7/b, n being an integer. 
By substituting 

Ex(r, y) = Eo(r) sin kyy, -Ho(r, y) = Ho(r) cos kyy, 


_in (3) and simplifying, equations in Es and Hj as functions of r alone are obtained. 


HH 
+ A(r) + B(r)He + eC(r)E, = 0, (6a) 
or? or 
+ A(r) —— + B(r)Es + uC(r) Hs = 0, (6b) 
or? or 
where the coefficients A(r), B(r) and C(r) have the values 
t(r) 4 1 22? 
2jRywr 
C(r) = J Rye J Rye 
Gr? (S2 + 
Y Equations (6a) and (6b) show that E» 


and HH are not independent in the bent 
pipe. Furthermore, E» and He do not 
vanish in this case, hence the methods 
of solution used for the straight pipe 
fail. Eg and Hy are not expressible in 
terms of Bessel functions. One possible 
method of solution of these equations, 
namely to substitute 


= (u)'/*He — 


and thus to obtain separated equations 
in ¢: and @2, is incorrect because the 
boundary conditions (4) are not satis- 
fied. 

Fic. 1 Equations (6a) and (6b) have been 
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solved completely by a method of approximation, using the theory of the 
Schrédinger equation with perturbations. Only the zero order and first order terms 
are considered. This does not affect the generality of the solution, because in 
practice the radius of curvature of the pipe, R, may be chosen very large compared to 
the constants of the equations and to the dimension a of the pipe. 

To rewrite (6a) and (6b) in the familiar Schrédinger form, let. 


(rHo)/R = B, 
r=R+s= R(1+5/R), 0<s<a, (8) 
y =2=/R. 
Thus 
0°B/ds? + fi(s)dB/ds + g(s)B + h(s)eEs = 0, (9a) 
0°Eo/ds? + fo(s)OEs/ds + g(s) Eo + h(s)uHe = 0. (9b) 


The coefficients, to the first approximation in R-, are given by 
fi(s) = R-(— 1+ 27°/K*), fa(s) = R-“(1 + 2y7/K?), 
g(s) = K* — ki — 2y°s/R = — 2y's/R, 


h(s) = 2jkyoy/RK?, 
where 


K? = 72+ = B+ 


Continuing the approximation, E, and 6 may be written as 


= (Eo + RUE) (10) 
and the perturbation of the angular coefficient for each case as 
R= + Re, for Ep, k2 = k?+ Rh, for B. (11) 


By substituting (10) and (11) in (9), the zero order and first order approximations 
may be written separately. For Eg these are 


0°(Es)o 


+ ki(Es)o = 0, (12a) 
O(E 
( (Es) o [1 + + (Esa 
Os? Os 
2j kyywou(He) 
+ €:(Es)o — 27%s(Ea)o + (12b) 


The zero order equation (12a) has the same form as the equation for the straight 
pipe, with the solution 
(Es)o = sin k,s, (13) 


where m and n are integers. 
Similar equations may be written for 8, giving 


r 
Bo = lim (He)o = Hmn COS kes = (He)o. (14) 
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Eq. (12b) may be rewritten as 


0*( 
ds K? 


This is the general form of the Schrédinger equation with perturbations, where the 
usual perturbation factor \ is equal to 1/R. 

By using the orthogonality condition for the Schrédinger theory, the value of e 
may be readily determined: 


f (— + 27's) sin? k,sds -f (1 + cos ks sin k,sds 
0 0 


f 2jkyopy mnEmn COS ks sin k,sds = 0. 
0 


Therefore 
= 


By using this value and (13), the first approximation (15) may be solved for (£):. 


The solution, satisfying the boundary conditions, is given by 


(Es)1 = Emn COS kas [(y*s)(a s)(2k.)*] 


+ sin kes {7 E + at | 


— Hm,» sin k,s[(jkysopy) J. (16a) 
In like manner, from the 8 approximation equations, 
hy = 


Since e;=/;, there is no change in the angle variable during the perturbations. The 
solution of the 8 equation, satisfying the boundary conditions, and corresponding to 
(16a) is 


B: = COS kes [(jkyswey)(keK?)] — sin ys [(jkywey) (BK?) ] 
+ Cos kes[(s/2)(1 — 2y2/K? + 7°/ki)] 
+ Hymn sin yas — 1 + y*s? + 2y2/K? — (16b) 


The complete solutions of (6a) and (6b), including both the zero order and first 
order approximations, may be written as 
Ey = {Emm Sin kys[1 + c1s] + Emn cos ks[s(a — s)ce] 
— sin ks[ucss]} {sin kyyei#*-20}, (17a) 
Hy = cos + cis] — Hm,» sin kys[s(a — s)co — 
— sin kes [ecs/ke] + Emn Cos k,s[ecss]} {cos kyye**'-2}, (17b) 
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where 
1 
2RR 2R RK? 2Rk, 
jhywy 7 1 
Rk, K? 2Rk, RK?k, 


By using (2) and the approximations (8) and (10), me components H,, H,, E, 
and £, are seen to become 


— K°H, = sin yk.(1 + cs — s/R + 
+ co(2s — a)y — jopecskys | 
+ COS kes[y{e1 + — cos(a — s)k,} 
+ sin kes[jweky(1 + cis + 2y%s/RK*) — kycsyes | 
+ Cos kes|jwekys(a — s)c2|}cos (18a) 
K*H, = sin — kes(a — + — 
+ Emn cos kes[jwek.(1 + cis + 2y2s/RK*) — jwe(2s — a)co + Rycsyes] 
— Hn,» sin ks [jopecs + yk,{s(a — — c4}] 
+ Hmn cos jwpekscss + yky(1 + cis — s/R + 2y%s/RK?) ]} 
‘sin (18b) 
— K°E, = sin — cos(a — s)ke + — 
+ Emn cos + cis — s/R + 2y2s/RK?) + jopekycss — cry(2s — a) ] 
— An,» sin + jouky{s(a — — cx} ] 


+ cos kes + cis + 2y2s/RK2) — 
‘sin (18c) 
— K°E, = { En,» sin kes[yky(1 + cis — s/R + 2y%s/RK*) — jones] 


+ Emn cos kes[ykys(a — s)co] 

+ sin kes[— jopk.(1 + crs + 2y*s/RK*) + — a) — kycssyu] 

+ Hmn cos kes[jou{c: — kes(a — + + jopkca)} 


> The solutions for the field components (18) satisfy the Maxwell field equations (1) 
within the approximation conditions imposed on the solution of the problem. 
For the special cases of Hn,, and En,, when one of the integers m or n is zero, 
the components may be obtained from (18). For m=0 and n not equal to zero: 


Es = Hon{jowykys(a — s)/RK»|sin kyye**™, 

He = Hon[1 — s/R — y*as?/2R + 7*s*/3R]cos kyye***™, 
= Ho,,[ys(a — s)R- + 7%) ]cos kyye®™, 
= Hon[RyyR7 (— wyea/K2 + R — y2as?/2 + y*s5/3) |sin 
—KnE, = Hon[joukyR + R — s — yas*/2 + y*s*/3)]sin 
E, = 0, 
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where 
Km = (K?) mao = ky + 

For m not equal to zero, n=0: 

E, = 0, 

Hy = {Hm,o cos kes[1 + cis] — Hm,o sin kes[(a — s)c2s — c4]} cos 
= { Hm,0 sin k,s[cxy(2s — a) — yk.(1 + as — s/R+ | 

+ cos + — cos(a — s)he} + kyye*™, 

H, = 0, E, = 0, ° 


| 
= 


— K2E, = {Hm,o sin kes[— jouk.(1 + cis + 2y’s/RK?2) + jopce(2s — a)] 
+ Hm,.o cos kes[jou{cr — kes(a — s)co + R} + |cos 
where the ¢1, C2, ¢3, C4 are calculated for m vanishing, and 
Ki, = (K*)nno = + 


It should be noted that both the Eo,, and En,o are missing. 

A consideration of the continued propagation of E and H waves from a straight 
pipe into a bent pipe yields some interesting results. A pure E»,, or Hn,, wave in the 
straight pipe will be reflected, partially, at the junction with the bent pipe. After 
reflection the amplitudes are proportional to a/R, and intensities to a?/R?, hence, for 
the first approximation, the reflected portion may be neglected. Thus a pure En, or a 
pure H,,,, wave in the straight pipe may be traced into the bent pipe, where it will 
become a mixed E and H wave. 

For a mixed £ and H wave in the straight portion, the intensities are proportional 
to a/R and must be considered. A mixed E and H wave in the straight pipe, because 
of the reflected portion at the junction, sets up an undetermined condition within the 
pipe, not predictable from the results of this paper. 

If the propagation constant is measured along the center line, a/2, of the bent 
pipe, there is no change in its value from that of the straight pipe. 
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AN APPLICATION OF ORTHOGONAL MOMENTS TO 
PROBLEMS IN STATICALLY INDETERMINATE 
STRUCTURES* 


BY 
W. M. KINCAID ano V. MORKOVIN 


Brown University 


1. Numerous methods have been devised for determining moment profiles’ in 
statically indeterminate structures. Most of these can be classified as either approxi- 
mate or exact. The approximate methods are usually simple to apply, but when em- 
ploying them it is necessary to specify numerical values for the dimensions and stiff- 
nesses (or their ratios) of the structure involved. The exact methods consist in solving 
systems of linear equations obtained by setting up relations between generalized dis- 
placements or by the equivalent means of using Castigliano’s Principle. As the degree 
of indeterminacy of the structure increases, the solution of the equations becomes 
more laborious. Most of the methods aim at reducing such labor by a suitable choice 
of unknowns. This paper is an attempt in that direction. 

2. Consider a structure whose degree of statical indeterminacy is N. Denote by 
M the true moment profile in the structure under a given load, and by My the moment 
profile under the same load when the structure has been made statically determinate 
by removing JN constraints (the so-called basic structure). The effect of the removed 
constraints may be replaced by the combined effect of N unknown generalized forces 
(couples and forces) Xi, X2, -- +, Xy. The generalized displacements of the loaded 
structure at the points of application and in the directions of X1, X2, ---, Xw, are 
assumed to be known and will be denoted by 61, 62, - - - , dy. (If X; is a couple, 4; is 
a rotation; if X; is an ordinary force, 6; is an ordinary displacement.) Let M; be the 
moment profile obtained when the force corresponding to X¥;=1, X;=0 for j #1, acts 
on the (unloaded) basic structure. Then the moment profile X;M; represents the effect 
of the ith constraint. Superposing the moment profiles due to the load and the con- 
straints yields the true moment profile M: 


N 
M = >, XiMi. (1) 
t=1 
The unknown quantities X; in (1) are to be determined by means of Castigliano’s 


Principle. 
Disregarding as usual the contributions of shearing stresses and axial forces, we 


get for the total strain energy U the following expression: 
1 
y= — (2) 


* Received July 1, 1943. 
! By the moment profile of a structure under a given load we understand the magnitude of the 
bending moment as a function of position; the graph of this function is commonly called the bending 


moment diagram. 
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where dx’ =dx/EI, E is Young’s modulus, J is the moment of inertia of any cross- 
section about the neutral axis, and the integral is taken over the entire structure. 
Castigliano’s Principle states that 


aU 
— = 45; (¢=1,2,---,W), 
or, by virtue of (1) and (2), 


N 
M.M + >> x, f = 6; (i =1,2,---,N). (3) 
j=t 

It will be observed that (1) and (3) are still valid if each X; represents a set of 
generalized forces (rather than a single force) acting simultaneously at different 
points. (See for instance Fig. 2d.) In such a case, each 6; would be made up of the 
generalized displacements corresponding to the given set of generalized forces. 

Previous attacks upon the problem have essentially consisted in choosing the 
points of application and lines of action of the unknown forces X; so as to make the 
system of equations (3) as simple as possible. Thus the moment profiles M; were 
completely determined. We propose to reverse this procedure by specifying the mo- 
ment profiles M; first. Let us choose these profiles so that 


f = 0 (i,j (4) 


Then each of the equations (3) will contain only one unknown, and we get at once 


6; — [MoM idx’ 
X;= (i= a°** (5) 
i 


The true moment profile M may now be obtained by substituting these values X; 
into (1). 

It will be recognized that equations (4) require that the moment profiles M; form 
an orthogonal system over the structure. Such a system can 
always be constructed by the standard orthogonalization 
process from the original set of moment profiles (or any simi- 
lar set of linearly independent moment profiles).? Therefore, 
the system of orthogonalized moment profiles and the corre- 
sponding generalized forces which appear in (5) will consist 
of linear combinations of the original sets of M; and Xj, re- 
spectively. In most practical applications, the exact form of 
these relations, as well as the physical interpretation of the 
generalized forces X;, is immaterial; only the orthogonalized 
moment profiles M; are needed. The simple example that 
follows will illustrate the notions introduced in the preceding Fic. 1 
discussion. 

3. Let us consider the triply indeterminate bent ABCD (Fig. 1), which has un- 
dergone a vertical displacement d and a rotation 8 at D (as alresult of a settlement of 


2 A discussion of orthogonal systems, explaining this process, will be found in R. Courant and D. 
Hilbert, Methoden der Mathematischen Physik (Berlin, 1931), vol. 1, p. 40, or in E. T. Whittaker and 
G. N. Watson, A Course of Modern Analysis (Cambridge, 1927), p. 224. 
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the foundation). The member 
UT UI AB is subjected to a uniform 


horizontal pressure p. The stiff- 
(2b) nesses EJ of the membersAB and 
| 


(2a) CD are assumed to be equal. In 
this and succeeding examples, 
the positive direction will be 
taken as downward in vertical 
members and toward the right 


in horizontal members. 
ines The solution will be carried 
out as follows: (1) the construc- 


tion of an orthogonal system will 
be shown, and (2) this system 
| 


Yh Yh 


(2d) will be used to obtain the desired 
moment profiles. 

Figures 2a, 2b, and 2c show 

STA iy, three self-equilibrating systems 
of generalized forces and the 

corresponding moment profiles; 
it is obvious that any set of re- 
actions (in equilibrium) can be 


built up by superposing these 

three in the right proportions. 
y (2e) 7 We note that the moment 
profiles Mz, and Mz are sym- 
K+2. A Kee metric with respect to a vertical 
i... = a axis through the midpoint of BC 
me ‘el and would be orthogonal to any 
antisymmetric profile. Therefore 


we replace M2. by Moa (Fig. 2d). 
It remains to replace M2, by a 
linear combination of Ma, and My, say aMo.+8Mw, that will be orthogonal to Mz. 
The orthogonality condition reads 


0= MM + ef = all’ + h’) + + 2H’), (6) 


where 1’ =1/EI gc, h’=h/EI4z. Eq. (6) will be satisfied if a=x+2 and B= —(x+1), 
where x=/’/h’. The resulting moment profile M2, is shown in 


Fic. 2 


v : Fig. 2e. The moment profiles Ma, Mea, and M2, form an orthogonal 
am system over the structure. From the manner in which this system 
: was derived, we see that it is independent of the loading and the 
Me choice of a basic structure. This is true not only with respect to 


this particular structure but in general. 
We choose a convenient basic structure and find the moment 


5 profile My due to the load (see Fig. 3). Eq. (1) now takes the 


phys form 


Fic. 3. M = Mo + + + (7) 
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Clearly, the quantities X may best be interpreted as mere coefficients by which the 
generalized forces in Figs. 2b, 2d, and 2e must be multiplied in order to yield, upon 
superposition, the correct reactions at the supports.* As for the interpretation of 
5a, dea, and 52, we recall that these quantities are equal to the partial derivatives, 
with respect to the corresponding quantities X, of the work done by the reactions. 
For instance, 


bea = 


2x 2d 
| Xeon + (« +1)X2.}0 — = 


2d 
Similarly —@ and 


Next, we evaluate the following important quantities: 


f = h'(x + 2), f Moadx! = +6), 


f Moedx’ = th'(x + 2)(2x + 1); (8) 
1 
f = f MoM 2adx' ff = + 2). (9) 


Thus the equations (5) specialize to 


— 66+ ph*h’ 60 — 12d/1 + 
6h' (x + 2) 2h'(x + 6) 
+ 10 — (3x + 2)p 


Sh’ (x + 2)(2x + 1) 


(10) 


To complete the solution we have only to substitute (10) into (7) and tabulate the 
values of M at A, B, C, and D. 

4. In the present method, the problem of solving the system (3) is replaced by 
that of constructing a set of orthogonal moment profiles. Once such a set is known for 
a given structure, the solution M corresponding to any loading is obtained merely 
by evaluating the right-hand member of (5) and substituting 
into (1).4 This is an advantage of the present method over others 8B Cc 
in which the complete system (3) has to be solved anew every 
time the load is altered. Thus, for important structures, it may 
be worthwhile to construct an orthogonal set, even in cases 
where the orthogonalization itself is fairly complicated. 

Fortunately, the rather lengthy standard orthogonalization A D 
process (see footnote 2) seldom needs to be used in its entirety, 
as the preceding example indicates. Shortcuts involving the use Fic 4. 
of symmetry, antisymmetry, and other characteristics of the 


’In most problems displacements at the constraints are zero. In such cases it is altogether un- 
necessary to visualize what particular generalized forces generate the orthogonal moment profiles because 
that physical notion is needed only for evaluating the 6;’s. 

4 The integrals JMjdx' which appear in the denominator of (5) are also independent of the loading 
and therefore can be computed once for all. 
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structure are usually available. For instance, in the case of the triply indeterminate 
frame in Fig. 4 (with I4g=Icp and Igc=I4p), we obtain the M;’s immediately from 
symmetry considerations (Fig. 5). 


Fie. 5. 


The evaluation of integrals [M Midx’ in (5) presents few difficulties, because 
the moments M; always vary linearly over any member and M, can be made to vanish 
over all but a few members by a suitable choice of the basic structure. Furthermore, 
there exist tables of such integrals for the common forms of moment profiles (trape- 
zoidal, parabolic, etc.).5 

Similarly, the remaining operations can be systematized so that even a person 
with little mathematical training can perform them. This is particularly true when 
numerical values of the stiffnesses EJ are known; then substituting (5) into (1) re- 
duces to taking scalar products on a computing machine. 

5. We conclude with an example illustrating an efficient arrangement of the work. 

Consider the four-legged bent shown in Fig. 6; assume that I4g=Icp=Jzr=Icu 
and Ipc=Icr=I re, and denote by the ratio 1’/h’ =lIaz/hI zc, as before. 


B Cc F G 
Me. 
A D E H 
Fic. 6. i Fic. 7. 


Since the structure is statically indeterminate of the ninth degree, we obtain nine 
orthogonal moment profiles Mi, M2, - - - , Ms, whose values at different points of the 
structure are given in the first ten columns of Table I. (We understand by M;(CB) 
the value of M; at the end of C of the member BC, and give a similar meaning to 
M;(AB), M;(BC), etc.) It will be observed that My, Me, M3, Mz, Ms, and Mg, are 
essentially reproductions of the moment profiles Ma, Mea, and M2, found for the 
two-legged bent in section 3. Considerations of symmetry are helpful in constructing 
M,, M;, and M,. The integrals {M?dx’ are now computed by means of the tables 
referred to in footnote 5, and are given in the eleventh column of our table. 

The parts of the table so far discussed can be used to find the moment profile due 
to any loading of the bent. Suppose the member BC is subjected to a set of vertical 


5 See, e.g., H. F. B. Miiller-Breslau, Die Graphische Statik der Baukonstruktionen (Leipzig, 1925), 
vol. 2, part 2, p. 56. 
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loads and moments. By selecting the basic structure as indicated in Fig. 7, we confine 
the moment profile Mo to BC. 
Defining 


1 1 5 
Py = Moxdx, = M,(l x)dx, 
we see at once that 


f MoM dx’ = I' [P,M<(BC) + (¢=1,2,---,9). (11) 


The values X; (twelfth column) are now obtained by dividing the negatives of the 
right members of (11) by the corresponding values in the eleventh column. 

To find the value of the true moment profile at (say) A, we multiply each term 
MAB) (first column) by the corresponding X; and add the results, and similarly 
for other points. It would ordinarily not be necessary to carry out the algebraic 
simplification of these sums, but this has been done for the sake of compactness: 


M(AB) = { (4x — 170x? — 414« — 225)P, — 2(« + 1)(10«? — 6x — 3)P2}, 

M(BC) = { — (64x + 628K? + 1188« + 585)P, + + 1)(16x? + + 51) Pi}, 
M(CB) = { (16x? + 11842 -+ 153x + 51) P, — + 1) + 14142 + 195« 72) Pr}, 
M(DC) = = { (20x? + 96x? + 233x + 135)P, + 2(— + Sta? + + 39) Pi}, 
M(CF) = {x(140x? 412x + 237)P, — 2(86x? + 3552 + 411x + 144) Pi}, 

M(FC) = = { — (124K? + 248x + 93)P, + + 179«? + + 36) Pi}, 
M(EF) = { x(36x? + 260x? + 377« + 135)P, — + 12542 + 176x + 69) Pi}, 
MGF) = { — + 208« + 135)P, + 2x(S4x? + 73x + 21) Pi}, 

M(FG) = { 2x(S4x? + 73x + 21)P, + 2(— — 38x? + Stx + 36)Pr}, 


1 
M(HG) = oT { «(36x® + 198k? + 284« + 135)P, — 2x(18x* + 90x? + 104% + 33)P,}. 


The arrangement given here is especially adapted to the use of computing ma- 
chines in case numerical values of x and the other quantities appearing are known.® 


6 After the manuscript of this paper had been completed, our attention was drawn to a work along 
similar lines by S. Miiller (S. Miiller, Zur Berechnung mehrfach statisch unbestimmter Tragwerke, Zentral- 
blatt der Bauverw. 1907, p. 23). Miiller introduces a “system of forces X;” instead of single forces, and 
reduces the system of equations to the diagonal form. However, his point of view and emphasis are quite 
different from ours, and his procedure is more lengthy than that presented here. 
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A STRAIN ENERGY DERIVATION OF THE TORSIONAL- 
FLEXURAL BUCKLING LOADS OF STRAIGHT COLUMNS 
OF THIN-WALLED OPEN SECTIONS* 


BY 
N. J. HOFF 
Polytechnic Institute of Brooklyn 


In the thin-walled open section columns of modern aluminum alloy aircraft tor- 
sional buckling and combinations of torsional and flexural buckling are of consider- 
able importance. The critical loads corresponding to these types of instability have 
been calculated by Wagner,! Kappus,? Lundquist and Fligg,? and Goodier* through 
integrating the differential equations of the problem. In the present paper the tor- 
sional-flexural buckling loads are determined with the aid of the Rayleigh-Ritz- 
Timoshenko method. This procedure obviates the derivation and integration of the 
differential equations as well as the geometric considerations connected with what 
Goodier termed “Wagner’s hypothesis.” 

The equilibrium of a straight bar of a length L and a cross-sectional area A, loaded 
axially with a compressive force of a magnitude oA distributed uniformly over the 
end section, can be investigated by assuming that each section of the bar undergoes 
a virtual displacement. The end sections of the bar are assumed to be restrained in a 
manner which precludes translations as well as rotations about any axis perpendicular 
to the end section, but which permits rotations about axes in the plane of the end 
section and warping of the end section. Barring displacements that would change 
the shape of the cross section (such displacements lead to plate- or shell-buckling), 
the most general virtual displacement pattern of the bar can be represented by the 
following infinite series: 


“= > a, sin (nwz/L), (1a) 
b, sin (nxz/L), (1b) 
B= > Cn Sin (nwz/L). (1c) 


* Received July 24, 1943. 

1 Wagner, Herbert, Verdrehung und Knickung von offenen Profilen, 25th Anniversary Volume of the 
Technische Hochschule, Danzig, 1929, p. 329. 
Torsion and buckling of open sections, N.A.C.A. Tech. Mem., No. 807, 1936. 

2 Kappus, Robert, Drillknicken sentrisch gedriickter Stibe mit offenem Profil im elastischen Bereich, 
Luftfahrtforschung, 14, 44 (1937). 
Twisting failure of centrally loaded open-section columns in the elastic range, N.A.C.A. Tech. 
Mem., No. 851, 1938. n 

? Lundquist, E. E., and Fligg, C. M., A theory for primary failure of straight centrally loaded columns, 
N.A.C.A. Report No. 582, 1937. 

* Goodier, J. N., The buckling of compressed bars by torsion and flexure, Cornell University Engineering 
Experiment Station, Bulletin No. 27, December, 1941. 
Torsional and flexural buckling of bars of thin-walled open section under compressive and bending 
loads, Transactions A.S.M.E., 64, A-103 (1942). 
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In these equations u is the virtual translation in the x-direction of a section of the 
bar at a distance z from the bottom section, v that in the y-direction, and 6 the virtual 
rotation of the section about its shear center. The x- and y-axes are parallel to the 
principal axes of inertia, and pass through the shear center of each section. The dis- 
placements of a section are shown in Fig. 1. 


Fic. 1. 


The increment 6U of the strain energy because of the virtual displacements is 


L L 
0 0 


L L 
+ 4GC f (d8/dz)%dz + f (d28/dz*) (2) 
0 0 


In Eq. (2) EI, and EJ, stand for the bending rigidities of the bar when the bending 
moment vector is parallel to the x-axis and the y-axis, respectively, and GC is the 
torsional rigidity calculated from the Saint-Venant theory of uniform torsion. The 
fourth term is due to the direct stress caused by non-uniform warping, and I is the 
warping constant defined in the theory of non-uniform torsion of thin-walled open 
sections. This theory is discussed in the previously mentioned references and in a 
paper written by the author.’ Eq. (2) follows from the theory of non-uniform torsion 
of thin-walled open sections and from the Bernoulli-Euler theory of bending, if the 
strain energy due to shear associated with bending and that associated with non- 
uniform warping are neglected. 


5 Hoff, N. J., Stresses in space-curved rings reinforcing the edges of cut-outs in monocoque fuselages, 
Journal Roy. Aeron. Soc., 47, 64 (1943). 
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Because of the orthogonal properties of the trigonometric functions the integrals 
indicated in Eq. (2) can be easily calculated: 


6U = + /4L*) >> 
n=1 n=1 


+ /4L*) ntc? + (x°GC/4L) >> nici. (3) 
n=1 n=1 

The decrease —5V of the potential of the external forces is equal to the work 6W 
done by them during the virtual displacements. The work déW done by the infini- 
tesimal force odA is equal to the force times the shortening of the distance between 
the end points of the fiber upon which it is acting. The shortening AL,,, of the dis- 
tance between the end points of the fiber passing through a point x, y can be calculated 
from the equation 


L 
AL zy af [(duz,y/dz)? +- (dv,,,/dz)? |dz. (4) 
0 
It may be seen from Fig. 1 that for small displacements 
Usy = u— JB, (5a) 
Vzy = 2 + xB. (5b) 


Upon substitution of the expressions of equations (1) and (5) into Eq. (4), integration 
yields 


AL = (?/4L) { naz + — 2y nance 
n=1 n=1 n=1 


+ nc + > meth (6) 


n=1 


The sum 6W of the work done by all the infinitesimal forces odA is 
= f o AL, ydA, (7) 
A 


where the integral is extended over the total cross-sectional area. With 
{sodA =P, the total compressive force, 
J/aydA =yoA, the static moment of the section with respect to the x-axis 
passing through the shear center, 
{axdA =x0A, the static moment of the section with respect to the y-axis 
passing through the shear center, 
Jfa(x?+y)dA = far’dA =I,, the polar moment of inertia of the section with respect 
to the shear center, and 
p?=I,/A 
equation (7) can be written as 


= (x*P/4L) { + > — 2y0 


n=1 n=1 n=1 


+ >> + > (8) 


n=1 n=1 
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According to the principle of virtual displacements the bar is in equilibrium in 
its original straight-line form if the change of the total potential 6(U+V)=d5U—6W 
is zero for any virtual displacement provided that first order small terms alone are" 
considered. Since both 6U (Eq. (3)) and 6W (Eq. (8)) contain only second order 
terms in the dp, by, Cn, the original straight-line form is a configuration of equilibrium. 
This equilibrium is stable only if the total potential increases, that is 6(U+V) is 
positive, for any virtual displacement. With the notation 


N = w*El,/L?, (9a) 
Q = wEI,/L?, (9b) 
R = (9c) 


the increment of the total potential can be written in the form 


6(U + V) = >> [n*(N/P) — 1]a2 + [n2(Q/P) — 


n=1 


+ [n*(R/P) + GC/P) — + 2yodncn — 2xobnen}. (10) 


With the notation 
A, = n*(N/P) 1, B, = n*(Q/P) 


C, = n?(R/P) + GC/P) — p’, (11) 
F, = — Xo, Gn = Vos 
X, = + + Cyc? + 2Fnbacn + (12) 


the infinite sum on the right hand side can be written as > \7_,X,. The necessary and 
sufficient conditions for its positive definite character are that all X, for m=1, 2, - - - 
must be positive definite. Necessary and sufficient conditions for this are 


n*(N/P — 1) >0, (14) 
n?(Q/P) — 1 > 0, (15) 
p?[n*(N/P) — 1][n*(Q/P) — 1] [(T/P) 1] 
— x3[n2(N/P) — 1] — yi[n*(Q/P) — 1] > 0. (16) 


Since inequalities (14)—(16) are necessary as well as sufficient conditions of stability, 
the bar may buckle if any one of them is not satisfied. Neutral equilibrium prevails, 
therefore, if any one of the following “buckling conditions” is fulfilled: 


n?(N/P) —1=0, (17) 
n?(0/P) —1=0, (18) 

p?[n*(N/P) — 1] [n*(Q/P) — 1][(7/P) 1] 
— x5[n*(N/P) — 1] — yo[n*(Q/P) — 1] = 0. (19) 


It is easy to prove that the original straight-line form of the bar corresponds to 
stable equilibrium if the compressive force P is sufficiently small, since by decreasing 
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P the first three terms on the right hand side of Eq. (12) can be made positive and 
large as compared to the last two terms. In the investigation of the stability of the 
bar under increasing values of P it is to be noted that if any of the inequalities (14)- 
(16) is satisfied for n=, it is also satisfied for n=k+ , where k and p are arbitrary 
positive integers. Consequently the bar is stable if, and only if, inequalities (14)—(16) 
are satisfied when »=1. The smallest critical load can be calculated from Eqs. (17)- 
(19) if 1 is substituted for m. This load alone is of practical importance unless geo- 
metric constraints (for instance rigid end fixation or additional supports between the 
ends of the bar) prevent displacement corresponding to the first terms of the Fourier 
series in Eqs. (1a)—(1c). In such a case the smallest value of that is compatible 
with the restraints must be used in Eqs. (17)—(19) for the calculation of the buckling 
load, for instance »=2 when the ends of the bar are prevented from rotating and 
warping. 

Eqs. (17)-(19) permit a discussion of the various types of buckling of bars of 
different cross section. With an asymmetric section x»#0 and yo#0. In this case 
with increasing P a value is reached at which Eq. (19) is fulfilled while the left sides 
of Eqs. (17) and (18) are still greater than zero. The buckling load P can be calculated 
from Eq. (19) which is a cubic in (1/P). The deflection pattern is flexural-torsional 
since it contains the non-vanishing coefficients @,, b,, and c, simultaneously. 

If the section has one plane of symmetry, one of the coordinates of the centroid, 


say Xo, vanishes. Then Eq. (19) reduces to 
[n*(Q/P) — 1] {p*[n*(W/P) — 1][(T/P)— 1] — = 0. (20) 


Consequently two distinct types of buckling are possible. One is purely flexural and 
symmetric. It corresponds to a buckling load which is the solution of Eq. (18). The 
other is flexural-torsional since it simultaneously contains displacement components 
corresponding to the non-vanishing coefficients a, and c,; it is antisymmetric; its 
buckling load P is the (smaller) root of the quadratic in (1/P) that can be obtained 
by dividing Eq. (20) by [n*(Q/P) —1]. The smaller of the two distinct buckling loads 
is of practical importance. The buckling load according to Eq. (17) is always greater 
than the smaller root of the quadratic. 

Finally, if the section is doubly symmetric or point symmetric, x»=yo=0. Then 
Eq. (19) reduces to 


(T/P) -1=0. (21) 


Buckling is either purely flexural or purely torsional. Of practical importance is the 
smallest of the solutions for P of the three Eqs. (17), (18), and (21). 
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—NOTES— 


USE OF SINE TRANSFORM FOR NON-SIMPLY 
SUPPORTED BEAMS* 


By A. G. STRANDHAGEN (Carnegie Institute of Technology) 


The problem of non-simply supported beams is approached by various mathemati- 
cal procedures. In certain applications several of the common methods are long and 
tedious. By employing the sine transform a certain ease can be claimed for most cases. 

The definition of the sine transform of a function y(x) in the interval (0, /) is 


S[y(x) ] y(x) sin (nwx/l)dx = o(n). (0<2<i;n=1,2,---) (1) 
0 

Recalling that the expression of a function y(x) in a Fourier sine series is’ 

y(x) = 7 b, sin nxx/I, (2) 
n=l 
where 
b, = f y(x) (3) 
0 


it becomes evident that the connection between the sine transform and the coefh- 
cients of the Fourier sine series is 


S[y(x)] = (4) 


Forms given by Eq. (2) and Eq. (3) are altered for the sake of convenience as follows: 


y(x) = o(n) sin (nwx/l), (5) 
where 
: 
o(n) = S[y(x)] -{ y(x) sin (nwx/I)dx. (6) 


For example, consider the sine transform of (d*y/dx?) in the interval (0, /); by 
definition 


S|d?y/dx?] = f sin (nwx/l)dx. (n= 1,2,---) 
Integrating formally by parts gives 


Likewise the sine transform of (d‘y/dx*) in (0, /) is: 


* Received August 9, 1943. 
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S[d*y/dx*] = — = [(— 


where v(n) in (7) and (8) is defined by Eq. (1). 
Consider a beam fixed at x=0 with axial loads P. The intensity of transverse 
loading is g(x), Fig. 1. The differential equation and boundary conditions are as 


follows: 
1. d4y/dx* + k*(d*y/dx*) = q(x)/EI, O<x<J 


2. y(0) = y(2) = 0, 
3. =0, = 9, 
g(x) =0 when O<2x< 
= 6(x) when b<x<e, 
= 0 when c<x<l, 


where 


and c>b. Let #=P/EI, and primes "| 

indicate differentiation with re- p “hhh p 
spect to x. Let S[y(x)]=v(n). 
Transforming d‘y/dx* and d*y/dx* —>| 

and g(x) and substituting (0) vy 

= y(l) =y’’(1) =0, there results 


+ (nx/l)40(n) — = f sin (nwx/l)dx. 


Solving tor v(m), where a?=(kl/7)?, 
+ 
n(n? — a?) n*(n* — a?) 


v(m) = — (1/x)*y’(0) sin (nwx/l)dx. (9) 
b 


Since y(x) =(2/)>.°_ v(m) sin nwx/l, then 
y(x) = — sin (nwx/I) 
n=1 — a’) 


sin 


+ sin (n#%a) (10) 
b 


a?) 
The remaining rea condition y’(0) =0 gives the following: 
1 
| (z’) si 11 


Since = (1/2a*)(1—a cot ma), then y’’(0) becomes 


(1 — wa cot ma) = —a 
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Further simplifications are possible in Eq. (12). Interchanging formally the integral 
and summation signi and summing, the following is obtained for y’’(0): 


(1/EI) {= k(l — x) 
sin kl l 


¥"(0) = 


(1 — racotra) J, 
where 


sin kl 


n(n? — a?) Qa? 
Thus by substitution of (13) in (10), 


mEI (1 — racot ra)J, 


2% =, sin(mrz/l) 
py n?(n? — a®) J A(x’) sin (nwx'/)dx’. (14) 


Knowing the variation of @(x) it is a matter of integration to obtain the required re- 
sults. Now suppose that P=0, i.e., the beam is under no axial loads, and subject to 
the same boundary conditions. Thus k =a =0 in equations (9), (10), and (11) and then 
3 


6(x’) sin 


y"(0) = 
Again interchanging formally the integral and summation sign, 


y = J, x’) a! (x x 2’, 


where 
Dd (1/n*) sin (nxx/l) = { 2(x/l) — .3(x/l)? + (x/l)*}. (0 < < 2.) 
n=l 

The equation for the elastic line becomes 


= — — 3(x/I)* + — I)(x! — 2)dx’ 


3 i 
n=l n* b 

To be sure, further summation in finite terms is possible, but this will lead to y(x) 
being defined in distinct intervals in (0, /), as in the solution furnished by the classical 
methods of differential equations; unquestionably, this is a disadvantage in engineer- 
ing computations. The above results, however, remain in the desired form, with one 
function y(x) in (0, /) regardless of the discontinuities of transverse loading. 

In like manner other boundary conditions may be imposed, and other beam prob- 
lems, such as beams on elastic foundations, can be solved. 


a 
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THE TREATMENT OF DISCONTINUITIES IN BEAM 
DEFLECTION PROBLEMS* 


By C. L. BROWN (Purdue University) 


The multitude of methods of determining the deflections of beams all stem from 
the fundamental differential equation 


M (1) 


where the abscissa x is measured along the axis of the beam, and y denotes the deflec- 
tion, M the bending moment, and EI the bending stiffness. The most obvious method 
of determining y is direct integration of (1). However, in most cases the right hand 
side of (1) is but sectionally analytical. A differential equation of the form (1) is then 
written for each section of the beam. When these equations are integrated, two con- 
stants of integration appear for each section. The evaluation of these constants of 
integration, though elementary, is extremely cumbersome. 

It is possible to avoid this sectionalizing treatment through the use of Heaviside’s 
unit step function, well known from operational calculus. This function is defined as 


follows: 
H.(2) for x <a, (2) 
ax) = 
i for #>6. 
It can readily be seen that 
H, if 
(3) 
HA, if a<b. 
Furthermore, if u(x) and v(x) are analytic, the continuous solution of 
d 
= u(x) + — a) 
dx 
is given by 
y= f u(t)d§ + Hao f + C, (4) 
0 0 
P 
where C denotes a constant of inte- 
gration. me ‘ 
The use of the unit step function = ® 
in the analysis of beams with concen- Mo 
trated loads is best shown by an ex- 
ample. The bending moment of the 5 ’ 
beam in Fig. 1 can be written as R, Re 


Fic. 1. 
d’y 
EI — = M = — M,+ — — a). 
dx? 


* Received Aug. 26, 1943. 
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Integrating according to (4) and taking account of the fact that slope and deflection 
vanish for x =0, we find 


dy 
EI = — Mix + 4Rix? — 4H,P(x — a)? (5) 
x 
and 
Ely = — + — 3H.P(x — a)*. (6) 


Both slope and deflection are zero at x=1. Thus, from (5) and (6) 
— Mil — + — = 0. 
Solving for M,; and R, 
M, = Pab?/I?, R, = Pb?(3a + 5) 


The deflection is then given by 


Pab? Pb*(3a+b P(x — a)’ 
2/? 
Pp In order to illustrate the use of the unit 
step function in the analysis of beams with sec- 
tionally constant moment of inertia, we con- 


sider the beam shown in Fig. 2. The reciprocal 


h cx of the bending stiffness can be written as 
a 
5 EI El, 
Fic. 2. and the bending moment as 


Px 


Thus, 
a? x 
El, [1+ — PH,(x 
dx? 2 


or, considering (3) 


fy Ps Ff 
EI, — = —+— [d.(x — a) + — PHi(x — 6) — PH.(x— 5), (7) 
dx? 2 2 
in which the second term on the right hand side has been written in this particular 
form in order to facilitate the application of (4). Integrating (7) we obtain 


ET Ps P E (x — a) 4+ oH (x 
6 2 


(x — 5)? (x — 6)? 


4 
| 
| 
3 M = — — PH,(x — 5). 
2 
it 
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The deflection is zero at x =0 and x=/. The application of these conditions to (8) 
yields 


P | 
Cy (6b — 3a*b + a’), C2 = 0. 


If, in addition to concentrated loads, the beam carries loads which are uniformly 
distributed over certain sections, the treatment is similar. 


A VARIATIONAL PRINCIPLE FOR A STATE OF COMBINED 
PLASTIC STRESS* 


By G. H. HANDELMANT (Brown University) 


In a recent paper! M. A. Sadowsky has stated a heuristic principle of maximum 4 
plastic resistance which he has applied to several states of combined plastic stress. 4 
The principle states that “among all statically possible stress distributions (satisfying 
all three equations of equilibrium, the condition of plasticity, and boundary condi- 
tions), the actual stress distribution in plastic flow requires a maximum value of the 
external effort necessary to maintain the flow.” W. Prager, in a contribution to the 
discussion of this paper’, has shown that the principle can be so interpreted as to lead 
to the correct differential equation for a beam under combined torsion and tension. 
This note is concerned with an accurate statement of the principle together with a 
proof of its validity for the case of a beam in a perfectly plastic state under combined 
torsion and bending by couples, the cross-section of the beam having an axis of sym- 
metry. Specifically, we shall prove the following variational principle for such a sys- 
tem. 

Among all statically possible stress distributions in a beam under a given torque 
(satisfying the equations of equilibrium, the condition of plasticity, and boundary 
conditions), the actual stress distribution when plastic flow occurs is the one for which 
the bending moment is stationary. 

Let us choose the coordinate axes in the following fashion. y lies along the axis of 
symmetry of the cross-section, z passes through the center of gravity of the cross- 
section and is parallel to the generators of the cylindrical beam, and x is perpendicular 
to y and z. We assume that the strain velocities, vz, vy, ¥:, are given by the same ex- 
pressions as in the case of an incompressible elastic material; i.e., 


v2 = — wyz + 
Vy = wxz — 30(x? — y? — 22), 


= ww(x,y) — Oyz. 


* Received Sept. 9, 1943. 

t This note was prepared at the suggestion of Professor W. Prager while the author was a participant 
in the Program of Advanced Instruction and Research in Mechanics at Brown University and was pre- 
sented to the American Mathematical Society on Sept. 12, 1943 under the title of On a principle of M. A. 
Sadowsky. 

1M. A. Sadowsky, Journal of Applied Mechanics 10, A-65 (1943). 

* Journal of Applied Mechanics 10, A-238 (1943). 
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w and @ are constants (w=angle of twist per unit length per unit time) and w(x, y) is 
an unknown function. The components of the strain velocity tensor, which are 


1 
’ 
2 Ox 


etc., will then be 


6 
€zz = €yy = 2 = by, 
w {dow w {Ow 


The components of the stress tensor will be denoted by o,2, Gzy, etc. The material 
is assumed to be incompressible and perfectly plastic. The stress-strain relations can 
therefore be written in the form 


(G22 — — (Ose — Cys: (1) 


In addition the following yield condition must be fulfilled: 
Tyy)? + + (22 + 6(02, + oi, + o2,) 6k?, (2) 


where & is a given constant. Eq. (1) will be satisfied if the stress tensor has the form 


ow 
Orz = = Ory = O, Ca * — 
Ox 
Ow 
= —— + I, — (3) 
dy 


p=p(x, y) being an unknown function of x and y. The equations of equilibrium for 


such a stress system reduce to 
OC zz Oo ys 


Ox oy 
Consequently we may introduce a stress function ku(x, y) such that 
dy Ox 


Since the surface of the beam is not stressed, «=const. on the surface. For con- 
venience, this constant may be taken as zero. Combining Eq. (4) and the yield con- 
dition, Eq. (2), we have 


where u,=0u/0x, uy=0u/dy. Therefore the function » must be given by 
—k 


Uz — ty) (5) 


‘ 
4 
4 
5 
| 
e 
i 
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By means of Eq. (3) and Eg. (4), we can compute 0dw/dx and 0w/dy in terms of the 
stress function u. This will yield 

ow k ow —k 


Ox wy dy wy 


Taking cross derivatives and subtracting one finds 


which becomes by virtue of Eq. (5) 
2w 
[yes(1 — — + ay — — — = 0. (6) 


We shall now determine the Euler equation for the variational principle stated 
previously. The bending moment and torque are given by 


Bending Moment = ff o2ydxdy = ff y(1 — uh — 


Torque = 2e ff udxdy. 


The domain of integration is the cross section of the beam. We note that the sym- 
metry assumption has been used in writing the bending moment in the form above. 
According to the usual procedure, let us form the function 


U2, Uy) = — ut — + ry, (7) 


where \ is an unknown constant. The Euler equation can then be written as 


a/a a/a 
Ox \Ouz Oy \Ouy Ou 


Substituting for g from Eq. (7), one finds 


[yue(1 — — + — [yuy(1 — — +d = 0. (8) 
Ox oy 


The unknown constant ) is determined by the fact that the torque is prescribed. On 
the other hand, w and @ can not be given arbitrarily but must be found from exactly 
the same condition. Consequently Eq. (8) is the same as Eq. (6), which proves the 
variational principle. 


3 
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ON THE ANTENNA PROBLEM* 
By S. A. SCHELKUNOFF (Bell Telephone Laboratories) 


The following remarks are made apropos to Brillouin’s recent discussion! of math- 
ematical difficulties involved in the retarded potential method of solving the antenna 
problem. The approximations involved in the actual solution of the final integral 
equation may be a source of far greater errors than the approximations in the equa- 
tion itself (such as 7A and 7B of Brillouin’s paper). Fig. 1 shows the first maximum 
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Fic. 1. The first maximum input resistance of center fed cylindrical antennas in free space as a func- 
tion of the average characteristic impedance: (1) is Hallén’s approximation to the solution of the integral 
equation for infinitely thin cylindrical shells, (2) is Schelkunoff’s approximation to the solution of Max- 
well’s equations for solid cylinders, (3) is Schelkunoff's approximation to the solution of Maxwell's equa- 
tions for infinitely thin cylindrical shells, (4) is Gray's recent approximation to the solution of the integral 
equation for infinitely thin cylindrical shells. Ke=120 (log 2//a—1) =602—120 ohms, where / is the 
length of each half of the antenna, Q is the parameter in Brillouin’s paper. The range of Ka in the figure 
corresponds to 38 <//a<1080 and 8.66<2< 15.33. 


* Received Nov. 20, 1943. 
1 Leon Brillouin, The Antenna Problem, Quarterly of Applied Mathematics, 1, 201 (1943). 
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resistance of center fed cylindrical antennas in free space. Curves (1), (3) and (4) 
refer to infinitely thin cylindrical shells (for which the antenna current vanishes at 
the ends); (1) is Hallén’s approximate solution of the integral equation, (3) is 
Schelkunoff’s approximate solution of Maxwell’s equations,? and (4) is Gray’s recent 
approximation to the solution of the integral equation.’ The difference between (1) 
and (4) is due solely to the difference in the methods of successive approximations. 
L. V. King’s approximation to the integral equation gives Rmax,1= 4000 for K,=420; 
this is considerably higher than even Hallén’s approximation. Curve (2) was calcu- 
lated by the same method as (3) but for solid cylinders. The “end” or “cap” capaci- 
tance was estimated as explained elsewhere;‘ this estimated capacitance is probably 
higher than the actual capacitance and in that region where (2) and (3) diverge the 
true curve is likely to be somewhat lower. In the case of hemispherical ends, curve (2) 
is raised still further. In all these curves it is assumed that there is no excessive 
localized capacitance in the vicinity of the input terminals; the effect of such capaci- 
tance is to lower these curves. Incidentally it means that one should never assume a 
“point generator” except when the input terminals aré tapered to mere points. 


2S. A. Schelkunoff, Electromagnetic Waves, D. Van Nostrand Co., New York, 1943, Chapter 11. 

3 Marion C. Gray, A modification of Hallén's solution of the antenna problem, Journal of Applied 
Physics, 15, No. 1, Jan. 1944. 

‘ Reference 2, page 465. 
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